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line.  230 


f 


Page 


Figure  C-3  Microfiche  layout  with  each  logarithm  of  flux 

(counts/sec)  as  ordinate  for  each  channel  and 
abscissa  the  UT  (sec)  every  180  sec  for  900 
sec  (16  samples/sec),  ALT  [km],  LONG  [deg-E], 
LAT  [deg],  B  [Gauss],  L  [  ],  ILAT  [deg],  MLT 
[hr].  E.G.  Eight  panels  (15  min  each)  of  12, 


33,  70  electrons,  8  panels  110,  160,  235  keV 

electrons,  etc.  P  =  Protons,  A  =  Alphas, 

P0  *  Proton  Detector,  BG  *  Background. 
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Figure  D-l 

CIRA  Reference  Atmosphere  of  depth  versus 

altitude  (1965). 
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Figure  D-2 

B  (z,E)  versus  altitude, 
r 

276 

xvlll 


ACKNOWLEDGEMENTS 


I  am  grateful  for  having  had  the  rare  opportunity  to  work  under 
the  supervision  of  someone  with  the  insight  of  Dr.  Richard  M.  Thorne. 

I  thank  Dr.  S.  V.  Venkateswaran  for  three  and  one-half  years  of  guid¬ 
ance,  inspiration,  and  friendship.  Dr.  George  L.  Siscoe  has  added 
significant  breadth  to  my  knowledge  of  the  space  sciences.  Conver¬ 
sations  with  Drs.  Charles  F.  Kennel  and  Margaret  Kivelson  were  helpful 
and  inspiring.  I  owe  special  thanks  to  Drs.  Thorne  and  Kivelson  for 
editing  the  manuscript  and  suggesting  style  changes. 

Finally  I  thank  Miss  Sherry  Lovell  for  typing  this  dissertation 
and  Mrs.  Dorothy  Deutsch,  Mr.  Ted  Doty,  Mrs.  Kathryn  Drake  and  Mr.  Wes 
Radlein  for  support  and  friendship.  While  a  graduate  student,  I  have 
been  supported  by  the  U.  S.  Air  Force  Institute  of  Technology. 

The  S3-2  and  S3-3  data  bases  were  provided  by  the  prime  investigator. 
Dr.  A.L.  Vampola,  The  Aerospace  Corporation,  ElSegundo,  Ca.  I  am  most 
grateful  for  many  helpful  conversations  with  him  and  Dr.  J.  Fennell  and 
for  their  hospitality  while  I  analyzed  the  data  sets. 


xix 


VITA 


November  n 
1 967 — E . S . , 
1  967--Com,;n 
1970-1971-- 
1 979--Ma jor 


,  1 945--Born,  Philadelphia,  Pennsylvania 
Saint  Joseph's  College,  Philadelphia 
ssioned  USAF 

K.S.,  University  of  Michigan 
,  USAF 


xx 


ABSTRACT  OF  THE  DISSERTATION 


Relativistic  Electron  Precipitation: 

An  Observational  Study 

by 
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Doctor  of  Philosophy  in  Atmospheric  Sciences 
University  of  California,  Los  Angeles,  1980 
Professor  Richard  M.  Thorne,  Chair 

The  object  of  this  dissertation  is  to  analyze  satellite  data 
to  investigate  the  morphology  of  strong  diffusion  relativistic  electron 
precipitation  (REP)  events  with  the  ultimate  purpose  of  gaining  insights 
into  the  precipitation  mechanism  and  determining  atmospheric  con¬ 
sequences.  Analysis  of  fourteen  (14)  months  of  magnetic  electron 
spectrometer  and  proton  telescope  data  from  the  Air  Force  S3-3  satellite 
from  July  1976  to  September  1977  is  presented.  Analysis  techniques 
and  data  catalogues  are  presented  in  detail.  Three  distinct  classes 
of  REPs  were  discovered:  an  overwhelming  majority  (96.5%)  exhibited 
precipitation  in  a  large  range  of  energies  below  1  MeV  for  electrons 
with  concurrent  proton  precipitation;  the  smallest  class  (1.3%)  also 
exhibited  concurrent  electron  and  proton  precipitation  but  the  electron 


and  proton  precipitation  but  the  electron  precipitation  occurred  c.-'y 
at  energies  above  1  MeV;  a  third  class  (2.2%)  had  no  associated 
proton  precipitation  and  electrons  precipitated  only  at  energies  belov, 
a  few  hundred  keV.  In  addition,  all  REPs  have  a  10-20%  frequency  of 
occurrence  that  linearly  increases  with  magnetic  activity,  show 
seasonal  variation  with  equinoxial  maxima  and  solstice  minima,  have 
mostly  small  latitudinal  (-0.5L)  and  large  longitudinal  (1 000-1 0,000km) 
extents,  and  occur  most  frequency  (9  out  of  10)  in  the  local  night¬ 
time  sector. 

A  theoretical  approach  to  the  REP  mechanism  is  presented. 
Interaction  of  the  radiation  belt  electron  and  ring  current  protong 
populations  with  electrostatic  and  electromagnetic  ion  cyclotron  and 
whistler  mode  chorus  waves  can  account  for  the  three  classes  of 
precipitation,  respectively.  Examples  of  each  class  are  presented 
in  detail . 

Finally,  an  investigation  of  the  energy  deposition  rates  from 
these  events  and  their  ultimate  effect  on  the  ionized  and  neutral 
atmosphere  with  emphasis  on  ozone  modification  is  presented.  The 
local  hydroxyl  concentration  can  be  increased  up  to  a  factor  of  10 
and  the  local  ozone  concentration  can  be  decreased  by  10-30%  by  REP 
events. 
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Introduction 


HF  radio  waves  can  propagate  from  transmitter  to  receiver 
via  refractive  interaction  with  the  ionosphere  (see  Fig.  1).  This 
propagation  is  known  as  forward  scatter  or  skywave  propagation.  The 
received  signal  can  be  diminished  (in  power)  if  the  electron  density 
below  the  scattering  altitude  is  enhanced  in  any  way.  Signal 
absorption  occurs  during  the  day  due  to  the  natural  increase  of  free 
electrons  below  the  scattering  altitude  by  solar  ionizing  radiation 
(1000-1 100°A) .  Absorption  can  also  be  increased  by  precipitating 
energetic  particles  (e.g.,  solar  protons)  that  reach  altitudes  below 
the  scattering  altitude.  One  well  known  class  of  absorption  event  is 
the  Polar  Cap  Absorption  (PCA)  events  which  occur  over  the  polar 
region  following  the  outburst  of  energetic  protons  during  flare 
activity  on  the  sun.  A  lesser  known  but  more  frequent  class  of 
events  involves  the  precipitation  of  relativistic  electrons  and  these 
will  be  the  subject  of  this  dissertation. 

Relativistic  Electron  Precipitation  (REP)  events  were  first 
observed  as  anomalous  daytime  absorption  of  forward  scatter  (more 
correctly  -  ionospheric  refraction)  HF  signals  by  Bailey  and  Pomerantz 
(1965).  The  detection  system  used  by  Bailey  and  Pomerantz  consisted 
of  three  paths  each  in  Alaska  and  Antarctica  operating  near  23  MHZ  with 
skywave  scattering  points  over  the  invariant  latitude  (i\)  range 
58-77°  (L  *  3.6-20.0).  The  scattering  altitude  was  estimated  to  be 
-70  km.  Peak  signal  absorption  typically  occurred  between  62-66° 
(L  =  4. 5-6.0);  namely,  within  the  region  of  long  term  trapping  for 
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FIG  1. 


Illustration  of  one-hop  skywave  propagation  via  the 
Ionosphere. 


2 


radiation  belt  particles. 


Since  solar  protons  could  be  ruled  out  due  to  the  absence 

of  absorption  over  the  higher  latitude  polar  cap  region  (A>  72°; 

L  >  10),  absorption  was  attributed  to  a  locally  precipitated  source 

of  particles.  These  particles  dissipate  most  of  their  energy  at 

altitudes  determined  by  their  range  in  the  atmosphere  and  thus  by 

their  energy  (Aikin  and  Bauer,  1965).  As  shown  in  Fig.  2,  protons 

require  at  least  several  MeV  while  electrons  require  only  a  few 

hundred  keV  of  energy  to  penetrate  below  70  km.  This,  however, 

immediately  excludes  the  main  bulk  of  auroral  electrons  (-1-10  keV) 

as  the  precipitating  particles  causing  absorption  and  indeed  visible 

auroral  activity  (monitored  by  all  sky  cameras)  was  found  to  have  no 

correlation  with  dayside  HF  absorption  (Bailey  et  al . ,  1970).  Since 

protons  of  several  MeV  energy  are  tenuous  beyond  (L  -  3.5)  and  do 

not  exhibit  outer  zone  integral  flux  variations  of  the  order  estimated 

-2  -1  -1 

for  strong  ionospheric  absorption  (-40  cm  str  sec  for  protons), 
Bailey  and  Pomerantz  concluded  that  trapped  electrons  at  energies 
greater  than  400  keV  were  the  cause.  In  situ  satellite  measurements 
have  also  demonstrated  that  relativistic  electrons  trapped  from 
A-  55-60°  (  L=  3-4)  can  exhibit  order  of  magnitude  flux  decreases 
during  disturbed  times  (Frank,  1965).  Furthermore,  since  the  maximum 
observed  absorption  intensities  occur  at  higher  L  values  than  the 
normally  expected  peak  intensity  of  trapped  electrons,  Bailey  and 
Pomerantz  (1965)  concluded  the  mechanism  operates  more  efficiently 
near  the  trapping  boundary. 
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Since  Bailey  and  Pomerantz  only  observed  the  events  during 
the  daytime,  they  assumed,  inconclusively,  that  the  precipitating 
mechanism  is  more  intense  on  the  dayside.  It  must  be  noted  that  this 
type  of  measurement  is  highly  frequency  dependent.  This  is  supported 
by  the  fact  that  as  the  forward  scatter  observing  frequency  changed 
from  50  MHz  in  the  early  1950's  to  35  MHz  pre-IQSY  (Int'l  Quiet  Sun 
Year)  and  finally  to  23  MHz  post  IQSY,  ability  to  detect  less  intense 
events  improved.  The  day/night  occurrence  probability  of  REPs  will  be 
further  discussed  as  we  review  subsequent  measurements. 

Other  types  of  observations  were  soon  employed  to  verify 
inferences  based  on  the  forward  scatter  results.  Riometers  which 
integrate  the  effect  of  precipitating  particles  (and  radiation)  of 
all  energies  were  used  to  confirm  signal  absorption  (Bailey  et  al . , 
1970).  All  sky  cameras  (observing  auroral  activity),  ionosondes 
(for  Es,  sporadic  E  layer,  detection),  magnetometers  (Bailey  et  al . , 
1966)  and  balloon  X-ray  measurements  (for  energy  spectrum  information) 
(Barcus  and  Rosenberg,  1966)  basically  confirmed  signal  absorption 
differences  between  daytime  REPs  and  nighttime  Es  (absorption  by 
precipitating  5-20  keV  electrons).  In  addition,  estimates  of 
precipitating  electron  spectra  for  REP  and  Es  conditions  for  quiet 
to  severely  disturbed  magnetic  conditions  were  constructed  from  X-ray 
and  forward  scatter  measurements  (Bailey  et  al . ,  1966;  Bailey  et  al . , 
1970).  Comparison  of  X-ray  and  forward-scatter-deduced  spectra  for 
the  same  event  revealed  the  marked  inability  of  the  forward  scatter 
system  alone  to  detect  the  presence  of  relativistic  electrons  due 
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FIG  2.  Atmospheric  penetration  depths  for  protons  and  electrons 
assuming  vertical  incidence  (adopted  from  Aikin  and 
Bauer,  1965). 
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to  the  large  uncertainty  in  scattering  altitude  (Rosenberg  et  al . . 
1972).  Spectra  derived  from  various  sources  will  be  compared  later. 

Another  drawback  to  estimating  REP  spectra  from  forward 
scatter  measurements  alone  is  the  uncertainty  of  absolute  absorption. 
The  signal  power  received  is  enhanced  by  increases  in  the  electron 
density  at  the  scattering  altitude  as  well  as  diminished  by  increases 
in  electron  density  along  the  signal  path  below  the  scattering 
altitude,  thus  an  accurate  measure  of  absorption  is  difficult. 

REPs  (and  PCAs)  appear  as  absorption  on  forward  scatter 
records  by  day.  At  night,  although  the  scattering  altitude  is 
higher  (by  10-20  km),  there  are  insufficient  free  electrons  below  the 
scattering  altitude  to  produce  signal  absorption  and  even  the 
strongest  precipitation  events  appear  as  signal  enhancements.  Rio- 
meters,  however,  continue  to  show  absorption  during  such  nighttime 
events  (Bailey,  1968). 

Two  Alaskan  paths  at  _A.=  65°  and  62°  (L  =  5.7  and  4.6)  were 
chosen  for  long  term  statistical  study  (Bailey,  1968).  For  the  period 
1  March  1964  to  31  December  1966,  REPs  producing  absorption  >  8  db 
occurred  on  319  days  with  maxima  at  equinoxes  and  minima  at  solstices. 
The  93  medium  and  large  REPs  (absorption  >  14  db)  for  the  same  period 
are  not  associated  with  sudden  commencement  (SC)  geomagnetic  storms, 
but  the  events  mostly  occurred  on  magnetically  disturbed  days,  lasted 
from  1-6  hours  (with  some  seeming  to  last  24  hrs)  and  were  mainly 
confined  to  the  12  hour  period  centered  on  local  noon. 
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After  1967,  the  system  detection  capability  became  erratic 


due  to  the  increasing  solar  activity.  The  system  was  terminated  for 
various  additional  reasons  early  in  1973  (Sullivan,  1979).  Attempts 
were  made  to  obtain  some  1969  forward  scatter  data  for  correlation 
and  calibration  with  OV 1-19  satellite  spectrometer  data  (see  Appendix 
A)  but  all  of  the  forward  scatter  data  except  for  the  first  five 
months  of  1969  had  been  lost  (Lincoln,  1978). 

In  1968-1971,  the  Norwegian  Defense  Research  Establishment 
(NDRE)  undertook  a  program  entitled  "Omega  Research  in  Norway" 

(Larsen,  1973).  This  program  monitored  the  VLF  (12.3  kHz)  short  hop 
(433  km)  path  between  the  transmitter  at  Aldra,^A_=  63.7°  (L  ~  5.1), 
and  the  receiver  at  Tromsfl,  j\.=  66.5°  (L  ~  6.3),  both  in  Norway,  for 
PCA,  substorm,  REP,  and  stratospheric  warming  research.  Phase 
recordings  of  the  VLF  signal  perpendicular  (to  propagation  plane)  and 
parallel  magnetic  components  provide  sensitive  detection  of  ionization 
variations  below  the  scattering  altitude  of  about  60-65  km  daytime 
and  80-90  km  nighttime  (Larsen,  1973).  A  diurnal  effect  is  seen  as 
an  advance  shift  in  baseline  of  60-80  ysec  from  local  night  to  day. 

REPs  (and  PCA/solar  flare  events)  are  readily  observed  both  day  and 
night  as  phase  advances  beyond  the  normal  undisturbed  levels. 

Thorne  and  Larsen  (1976)  compared  the  Norwegian  VLF  and 
Alaskan  HF  observations  of  REPs  and  Auroral  Electrojet  (AE)  Indices 
(for  substorm  identification)  for  the  period  January-May,  1969.  They 
indicate  "...a  direct  correlation  between  the  onset  of  intense  night¬ 
time  REP  events  with  substorms  while  dayside  REPs  are  far  less  frequent. 
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are  shorter  in  duration  and  are  often  delayed  by  as  much  as  several 
hours  from  the  onset  of  substorm  activity"  (Thorne  and  Larsen,  1976). 
Figure  3  shows  for  example,  a  simultaneous  increase  in  AE  and  VLF  phase 
advance  in  Norway  at  2100  UT  followed  1-1  1/2  hrs  later  by  the  HF 
absorption  in  Alaska.  This  sequence  agrees  with  the  observation  of 
the  25-26  August  1967  substorm  by  Hones  et  al . ,  (1971)  but  not  their 
inference  of  gradient  drifting  electrons  injected  at  substorm  onset 
since  the  gradient  drift  lifetimes  of  relativistic  electrons  are  on 
the  order  of  minutes  (Thorne  and  Larsen,  1976).  In  addition,  only  a 
select  class  of  substorms  can  trigger  dayside  REPs  (Rosenberg  et  al . , 
1972;  Thorne  and  Larsen,  1976)  and  substorm  intensity  is  only  weakly 
correlated  in  the  selection  of  which  substorm  will  spawn  a  REP 
(Thorne  and  Larsen,  1976). 

Many  observers  using  riometers  (Brown,  1964;  Jelly  and  Brice. 
1967)  and  chains  of  riometers  (Ansari ,  1965;  Lukkari  et  al . ,  1975; 
Lukkari  and  Kangas,  1976;  Lukkari  et  al . ,  1977)  have  observed  similar 
relations  between  subauroral  energetic  electron  precipitation  and 
slowly  varying  absorption  (SVA),  substorm  onsets,  and  magnetic 
pulsations.  Riometers  observe  total  absorption  occurring  both  above 
and  below  the  scattering  altitudes  of  VLF  and  HF  systems  and  alone 
cannot  determine  where  (in  altitude)  the  absorption  occurs  nor  yield 
a  precise  precipitation  spectrum  (Bailey,  1968). 

Balloon-borne  measurements  of  X-rays  during  REP  events  (usually 
with  uncollimated  Nal  detectors  near  32  km  altitude)  have  the  advantage 
over  riometer  measurements  in  sensitivity,  time  resolution,  and 
the  fact  that  they  provide  an 
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FIG  3.  Time  sequence  of  (A)  HF  absorption  (dark  shaded  area) 
in  Alaska;  (b)  VLF  phase  advance  (shaded)  from  normal 
nighttime  trace  (dashed  line);  (c)  AE  index  for  13-14 
March  1969  from  Thorne  and  Larsen,  1976. 


energy  spectrum  (Bewersdorff  et  al.,  1968).  In  addition  to  the 
disadvantages  of  imprecise  balloon  location  (~  100  km  horizontally) 
and  proton  contamination  of  X-ray  scintillators,  the  theoretically 
predicted  X-ray  spectrum  (the  scintillator  measures  mostly  secondaries) 
is  very  sensitive  to  the  assumed  incident  electron  spectrum,  usually 
in  exponential  form  («*  exp  E/Eo),  and  atmospheric  depth  (Bailey  et  al . , 
1970).  Only  six  or  seven  X-ray  measurement  deduced  electron  energy 
spectra  during  energetic  electron  precipitation  events  have  published 
and  range  from  spectra  constructed  by  the  addition  of  soft  (Eo  =  9  KEV) 
and  hard  (Eo  =  40  KEV)  components  (Bailey  et  al . ,  1970)  to  very  hard 
(Eo  =  180-200  KEV)  spectra  (Rosenberg  et  al . ,  1972;  Parks  et  al . , 

1979). 

Rockets  have  carried  spectral  and  integral  electron  spectro¬ 
meters  (Matthews  and  Simons,  1973;  Vij  et  al . ,  1975;  Matthews  et  al . , 
1976)  as  well  as  having  boosted  parachute-borne  X-ray  scintillators 
(Vij  et  al.,  1975)  during  REP  and  auroral  events.  Simultaneous 
electron  (rocket  to  275  km  apogee)  and  X-ray  (rocket  boost  to  60  km) 
observations  show  fair  agreement  at  least  in  the  much  "softer" 
auroral  events  (Vij  et  al . ,  1975).  The  spectrometers  yield  a  more 
precise  electron  energy  spectrum  and,  depending  on  instrument 
orientation  and  flight  configuration,  add  pitch  angle  distribution 
information.  An  obvious  drawback  of  this  type  of  experiment  for 
REP  research  is  the  lack  of  morphological  statistics  due  to  prohibitive 
logistics  expense. 
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Satellite  X-ray  spectrometer  experiments  provide  simultaneous 
measurements  over  a  wide  range  of  local  times  because  of  their 
extensive  fields  of  view  (Imhof  et  al . ,  1974)  as  well  as  opportunities 
for  coordinated  studies  with  particle  spectrometers  for  improved 
spatial /temporal  resolution  (Imhof  et  al . ,  1975).  Comparison  of 
X-ray  and  particle  measurements  show  agreement  to  only  a  few  hundred 
keV  and  that  impulsive  (<  5  min  duration)  events  near  the  energetic 
trapping  boundary  were  substorm  correlated  (Imhof  et  al . ,  1978). 
Coordinated  two  satellites  X-ray  and  particle  measurements  show 
significant  variations  in  local  time  dependence  of  energetic  electron 
precipitation  though  fluxes  seemed  greater  in  late  morning  than  early 
morning. 

The  majority  of  satellite  REP  measurements  have  been  made 
with  integral  (Brown  and  Stone,  1972;  Larsen  and  Thomas,  1974; 

Morfill  et  al . ,  1975)  and  spectral  (Vampola,  1969;  Vampola,  1971 
a  band  b;  Imhof  et  al . ,  1973;  Imhof  et  al . ,  1975)  particle  spectro¬ 
meters.  REP  precipitation  is  observed  to  be  a  narrow  (AL  ~  .5-1.5) 
phenomenon.  Figure  4  shows  an  invariant  latitude  versus  magnetic 
local  time  plot  of  electron  precipitation  spikes  with  E  >  425  keV 
from  the  last  five  months  of  1967  from  the  low  polar  orbit  satellite 
0G0  4  (Brown  and  Stone,  1972).  Notice  a  relative  scarcity  of 
events  between  06-12  MLT.  Also,  the  events  extend  from  approximately 
the  middle  of  the  outer  electron  zone  to  the  limit  of  closed  field 
lines  with  the  highest  density  clustered  about  the  limit  of  stable 
trapping.  A  general  shift  to  lower  invariant  latitudes  is  observed 
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with  increasing  geomagnetic  activity.  Event  time  scales  range  from 
several  minutes  to  several  hours.  Approximately  five  times  as  many 
events  occur  at  2200  MLT  as  at  1000  MLT.  This  study  was  extended 
by  the  use  of  the  >  500  keV  measurements  on  board  another  low  polar 
orbit  satellite  ESRO  4  with  similar  results  (Worfill  et  al . ,  1975). 

Another  investigation  of  magnetic  electron  spectrometer 
data  of  precipitating  electrons  in  the  energy  range  225  keV  <  E  < 

2.47  MeV  aboard  polar  orbit  satellite  0V3-3  for  the  period  August 
1966-September  1967  show  common  occurrences  of  precipitation  at  local 
evening  and  midnight  in  a  narrow  region  (AL  ~  0.5)  at  the  local  limit 
of  trapping  (Vampola,  1971).  Rare  events  (3)  at  local  evening  just 
inside  the  plasmapause  where  higher  energy  electrons  (>  1  MeV)  are 
preferentially  precipitating  were  also  found. 

Measurements  of  precipitating  electrons  from  >  48  keV  to 
>  434  keV  on  the  polar  orbit  satellite  ESRO  1A  during  the  REP  of 
2  February  1969  at  local  morning  at  /V=  63°-70°  (L  =  4. 9-8. 5)  have 
been  reported  by  Larsen  and  Thomas  (1974).  The  energy  spectrum  for 
this  event  is  a  bit  "harder"  than  those  previously  obtained  by  rocket 
measurements  (Matthews  and  Simons,  1973)  but  much  harder  than  those 
deduced  from  balloon  X-ray  measurements  (Bailey  et  al . ,  1970). 

Figure  5  compares  REP  energy  spectra  from  various  sources 
previously  discussed  with  typical  discrete  and  diffuse  auroral  spectra 
(Lui  et  al ■ ,  1977).  Satellite-borne  spectrometers,  particularly 
spectral  spectrometers,  have  measured  the  hardest  spectra  in 
comparison  with  non  satellite  spectra  and  represent  the  most  direct 
means  of  examining  REPs  of  all  intensities. 
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FIG  5.  Differential  energy  spectra  of  precipitating  electrons  of 

the  following  types:  (A)  Diffuse  aurora  (Lui  et  al . ,  1977); 
(B)  Discrete  aurora  (Lui  et  al. ,  1977);  (C)  Balloon  X-ray 
(Bailey  et  al . ,  1970);  (D)  Rocket  (Matthews  and  Simons, 

1973);  (E)  Satellite  (Larsen  and  Thomas,  1974);  (F)  Satellite 
(Vampola,  1971). 
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It  was  the  purpose  of  this  dissertation  to  determine  the 
mechanisms  responsible  for  REP  events,  their  relative  importance 
and  atmospheric  consequences  through  the  study  of  satellite  data. 
Magnetic  electron  spectrometer  data  on  the  low  polar  orbit  satellites 
0V1-19  and  S3-3  were  graciously  provided  by  the  instrument  builder 
and  prime  investigator,  A.  L.  Vampola,  Aerospace  Corp.,  El  Segundo, 

Ca.  An  exhaustive  study  of  the  0V1-19  data  set  was  performed  but 
for  various  reasons  the  data  set  yielded  only  marginally  significant 
results.  This  study  is  described  in  detail  with  results  in  Appendix  A. 
The  S3-3  data  set  has  few  of  the  drawbacks  of  the  0V1-19  data  set  and 
is  the  main  source  of  study  here.  The  S3-3  satellite  orbit  and 
orientation  as  well  as  particle  and  wave  measurements  are  described 
in  Appendix  B. 

This  thesis  is  arranged  with  four  main  chapters.  In  Chapter  1, 
details  of  the  S3-3  REP  event  morphology  are  presented.  A  description 
of  the  objective  and  subjective  analysis  of  the  data  needed  to 
identify  events  is  contained  in  Appendix  C. 

Chapter  2  contains  a  general  background  of  the  types  of 
mechanisms  previously  postulated  and  others  newly  discovered  capable 
of  explaining  REPs.  The  mechanisms  are  compared  with  the  classes  of 
events  found  in  Chapter  1.  In  addition,  a  case  study  of  each  type 
is  presented  with  the  view  of  advancing  the  knowledge  of  each 
mechanism. 

Chapter  3  compares  the  atmospheric  energy  deposition  profiles 
for  REP  and  REP  X-rays  with  profiles  from  other  sources. 
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Profiles  for  REPs  of  various  classes  and  intensities  are  presented 
as  well  as  determinations  of  their  importance  and  consequences  in 
ionization  production  and  ozone  chemistry.  The  electron  and  X-ray 
deposition  calculations  are  presented  in  Appendix  D. 

Chapter  4  summarizes  the  findings  and  discusses  the 
conclusions  and  outlooks. 
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CHAPTER  1  .  EVENT  MORPHOLOGY 


1 .1  Results 

We  present  here  the  statistical  data  on  REP  occurrence  measured 
on  S3-3  during  a  14  month  period  from  11  July  1976  to  29  September  1977. 
The  S3-3  data  collection  specifications  are  described  in  Appendix  B. 

The  methodology  of  selecting  REPs  as  well  as  the  data  base  description 
is  provided  in  Appendix  C.  Fig.  1.1,  a-c,  shows  a  summary  of  the  S3-3 
data  set;  the  plot  includes  daily  values  of  the  highest  energy  channel 
(in  keV)  exhibiting  strong  diffusion  precipitation  (Panel  A);  the 
daily  number  of  outer  zone  (L=4-8)  passages  (Panel  B);  the  daily  mean 
of  the  hourly  equatorial  DST  in  gamma  units  (Panel  C);  and  the  daily 
Ap,  equivalent  daily  planetary  magnetic  amplitude  in  2  gamma  units 
(Panel  D).  Outer  zone  coverage  (Panel  B)  was  fairly  uniform  through¬ 
out  the  data  set  (although  29  days  had  no  coverage)  averaging  a  little 
over  5  outer  zone  passages  per  day.  Coverage  was  a  function  of  tape 
recorder  capacity,  power  limitations,  microfiche  availability,  etc. 

DST  and  Ap  are  presented  to  represent  the  relative  magnetic  condition 
of  the  magnetosphere  since  AE,  the  Auroral  Electrojet  Index,  was  not 
yet  available  for  this  time  period.  AE  indices  for  five  stations,  AE 
(5),  for  seven  2  day  periods  in  1977,  however,  were  available  and  are 
discussed  in  the  next  section.  (Allen,  private  communication,  1979.) 

1 .2  Frequency  of  Occurrence  and  Magnetic  Dependence 

The  monthly  average  frequency  of  occurrence  of  REP  events  ob¬ 
served  during  the  satellite  transit  of  the  outer  radiation  zone  (L=4-8) 
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FIGURE  1  .la, 


I 
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b,c.  Daily  values  from  11  July  1976  to  29  September  1977. 
Plots  are  given  of  REPs  with  highest  energy  channel 
exhibiting  strong  diffusion  (Panel  A);  the  total 
number  of  outer  zone  passings,  L=4-8  (Panel  B); 
daily  mean  of  the  hourly  Dst,  y  (Panel  C);  Ap, 
equivalent  daily  planetary  magnetic  index,  2y  (Panel 
D). 
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is  represented  in  Fig.  1.2.  Maximum  occurrence  is  found  near  both  (1976, 
1977)  autumnal  equinoxes,  at  and  si  ightly  beyond  vernal  equinox,  and  near 
winter  solstice  while  pronounced  minima  are  evident  near  both  (1976, 
1977)  summer  solstices.  This  is  generally  consistent  with  the  equi- 
noxial  maxima  and  solstice  minima  found  by  Bailey  (1968).  In  this 
sense,  the  winter  maxima  is  anomalous,  but  approximately  five  strong 
magnetically  disturbed  periods  took  place  in  Dec  1 976- Jan  1977  (see 
Fig.  1.1b)  probably  accounting  for  the  rise  in  activity.  The  magnetic 
dependence  is  apparent  when  the  same  monthly  percentage  of  REPs  nor¬ 
malized  to  outer  zone  coverage  is  plotted  for  days  when  Ap  >  30 
(Panel  A),  20  <_  Ap  <  30  (Panel  B),  10  <_  Ap  <  20  (Panel  C),  and  Ap  <  1 0 
(Panel  D)  in  Fig.  1  .3  The  anomalous  winter  maximum  in  Dec  is  apparent 
for  30  <  Ap  <  20  while  the  Jan  maximum  is  seen  only  for  20  <  Ap  <  30. 

The  remaining  maxima  and  minima  are  fairly  consistent  with  magnetic 
conditions. 

Plotting  the  REP  occurrence  frequencies  for  all  14  months  for 
different  ranges  of  Ap  (magnetospheric  disturbance)  (Fig.  1.4),  one 
sees  a  nearly  linear  increase  in  REP  occurrence  with  magnetic  distur¬ 
bances. 

Fig.  1.5  shows  S3- 3  coverage  of  seven  2  day  periods  by  AE  index 
(Allen,  private  communication,  1979).  This  data  is  merely  presented  to 
support  the  fact  that  not  all  substorms  trigger  REPs  (note  particularly 
Panels  A  and  E)  in  agreement  with  the  findings  of  Rosenberg  et  al . 

1972)  and  Thorne  and  Larsen  (1976).  The  local  times  of  collections  for 
both  Panels  A  (1  outer  zone  pass  per  rev)  and  E  (2  outer  zone  passes 
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Monthly  REP  event  percentage  of  occurrence  normalized 
to  outer  zone  coverage  for  July  1 976-September  1977. 
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FIGURE  1.3  Monthly  percentage  of  REPs  normalized  to  outer  zone 
crossings  (L=4-8)  for  magnetic  conditions:  Ap  >  30 
(Panel  A);  20  _<  Ap  <  30  (Panel  8);  1  0  _<  Ap  <  20  (Panel 
C);  Ap  <  10  (Panel  D).  Values  for  JAN  and  FEB  in  Panels 
A  and  B  overlap.  P=0  for  both  JAN  and  FEB  in  Panel  A 
and  P=0.86  and  0.39  for  JAN  and  FEB  respectively  in 
Panel  B. 
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FIGURE  1.4  REP  event  percentage  of  occurrence  normalized  to  outer 
zone  coverage  (L=4-8)  versus  magnetic  conditions  for  all 
months  JUL  1976-SEP  1977.  Total  number  of  outer  zone 
passes  for  each  category  are  in  parentheses. 
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FIGURE  1.5  Seven  2  day  periods  of  AE  index  (y)  throughout  1977 
with  REPs  indicated  by  arrows  (t)  and  hours  of  noncov 
erage  by  shaded  horizontal  bar  (■■).  AE  =  AU  (Upper 
Trace)  -  AL  (Lower  Trace). 
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per  rev)  occurred  near  midnight  where  one  would  expect  highest  prob¬ 
ability  of  detection  (see  Fig.  1.7  for  example).  Many  of  the  outer 
zones,  particularly  the  outer  edge,  in  the  passes  of  Panels  A  and  E 
were  very  irregular  showing  trapped  flux  variations  to  nearly  an  order 
of  magnitude  though  no  precipitation  was  evident  at  the  time  of  the 
pass.  Also  noticeable  are  times  when  substorms  with  large  values 
(~  900y )  of  AE  do  not  spawn  REPs  (03-04L,  30  Jan  77,  Panel  A;  12-14L, 

29  Jul  77,  Panel  E;  03-1 1L,  30  Jul  77,  Panel  E;  13-18L,  19  Sep  77, 

Panel  F)  while  substorms  corresponding  to  smaller  AE  may  be  associated 
with  REPs  (08-1 OL,  10  Mar  77,  Panel  B;  01-02L,  8  Apr  77,  Panel  C; 
05-06L,  3  May  77,  Panel  D;  04-06L,  20  Sep  77,  Panel  F;  11-12L,  23  Sep 
77,  Panel  G).  This  agrees  again  with  Thorne  and  Larsen's  (1976)  find¬ 
ing  of  weak  correlation  between  AE  index  and  REP  observations.  An 
exception  to  this  correlation  is  the  REP  activity  at  03-04L,  22  Sep  77, 
Panel  G,  where  a  1500  y  storm  has  associated  REP  activity  for  all  four 
outer  zone  passes  on  the  same  revolution  (see  event  catalogue.  Appendix 
C,  REV  P  3569).  This  orbit  was  oriented  in  the  0200  and  1400  hrs  MLT 
plane  and  could  infer  a  longitudinal  extent  of  nearly  10,000  km. 

1 .3  Local  Time  Dependence  and  latitudinal  Extent 

To  examine  the  local  time  morphology,  REPs  (Fig.  1.1,  Panel  A) 
were  plotted  for  a  given  date  versus  magnetic  local  time  (MLT)  in  Fig. 
1.6.  In  order  to  test  for  any  bias  in  the  excursion  of  the  orbit  plane 
with  time  (skewed  solid  lines),  we  have  sketched  (hatched  areas)  the 
span  of  local  times  over  which  REPs  could  be  observed  in  the  range 
4  <  L  <  8  for  a  given  local  time  orientation  of  this  orbit  plane.  A 
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FIGURE  1.6  REPs  by  date  versus  MLT  (hrs).  Hatched  area  indicates 
regions  of  coverage  (4  L  _<  8)  around  the  local  time 
of  the  orbit  plane  (skewed  solid  line). 
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possible  apogee-perigee  bias  in  the  orbit  plane  was  investigated  but 
after  careful  analysis  (see  Appendix  C)  the  bias  was  found  to  be  insig¬ 
nificant.  The  most  striking  feature  is  the  pronounced  night-day  asym¬ 
metry  with  the  greatest  number  of  REPs  occurring  in  the  night  sec. or. 

A  plot  of  the  ratio  of  day  to  night  REPs  is  presented  in  Fig.  1.7  as 
a  function  of  magnetically  disturbed  conditions.  Overall  there  is  a 
9:1  preference  for  nighttime  REP  events.  This  local  time  asymmetry 
is  less  for  highly  disturbed  (Ap  >  30)  conditions. 

The  actual  local  time  morphology  of  REP  events  is  shown  in 
Fig.  1.8  (dots  indicate  very  narrow  localized  events  and  lines  are 
drawn  for  extended  events)  where  they  are  plotted  against  invariant 
latitude  for  varying  magnetic  conditions.  Most  events  show  narrow 
latitudinal  extends  -0.5L  (dots)  while  a  few  showed  excursions  up  to 
4-5L,  mostly  in  the  235,  435  keV  electron  channels.  Readily  apparent 
is  the  relative  lack  of  dayside  REPs  particularly  during  magnetically 
quiet  days  (Ap  _<  10).  Furthermore,  the  scarcity  of  REPs  over  the  local 
time  range  1300-1900  remains  apparent  for  Ap  ^  10.  As  can  be  seen 
from  Fig  1.6  the  orientation  of  the  orbit  plane  for  these  MLTs  cov¬ 
ered  roughly  all  months  except  the  Mar-May  maxima.  Both  autumnal 
equinoxes  and  the  anomalous  Dec-Jan  maxima  were  covered  so  we  believe 
that  the  statistics  for  the  local  afternoon  hours  are  not  biased  by 
the  orientation  of  the  orbit  plane:  thus  the  tenuous  afternoon  statis¬ 
tics  appear  to  be  real.  Also,  for  highly  disturbed  conditions  ( Ap  >_  20) 
the  event  location  moves  towards  lower  invariant  latitudes;  this  is 
most  noticeable  for  Ap  >  30. 
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FIGURE  1  .7 


Ratio  of  day  (MLT  06-18)  to  night  (MLT  18-06)  for 
varying  magnetic  conditions. 
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FIGURE  1.8  REPs  displayed  by  invariant  latitude  (degrees)  versus 
MLT  (hrs)  for  varying  magnetic  conditions.  Lines  in¬ 
dicate  excursions  of  events  in  invariant  latitude. 
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1 .4  Longitudinal  Extent 


The  slow  (1.5  hr/month)  excursion  of  the  orbit  plane  in  local 
time  (see  Appendix  B)  does  not  allow  direct  analysis  of  longitudinal 
extent.  Estimates,  however,  can  be  made  from  previously  published 
ground  based  events  (e.g.  Bail ey,  1968)  though  these  would  probably  be 
a  lower  limit.  In  addition,  general  limits  could  be  estimated  on  an 
individual  case  basis  from  the  satellite  revisit  times  (~3  hrs)  where 
consecutive  coverage  exists. 

Duration  of  ground  based  REP  events  was  observed  to  be  1-6  hrs 
with  some  events  lasting  for  24  hrs  (Bail ey,  1968).  Several  S3-3 
REPs  seemed  to  last  for  24  hrs  out  continuity  of  these  events  is  un¬ 
certain  since  events  can  occur  over  a  period  of  many  substorms  (e.g. 
Fig  1.5,  Panel  F).  The  maximum  longitudinal  extent  at  60°  latitude 
is  approximately  20,000  km.  Bailey's  1  hr  lower  limit  translates  into 
850  km  which  roughly  agrees  with  Larsen's  (1973)  estimate  of  REPs  ex¬ 
tending  greater  than  1000  km.  The  6  hour  duration  estimate  translates 
to  a  5,000  km  longitudinal  extent.  For  very  disturbed  times,  all  four 
outer  zone  crossings  on  the  same  3  hour  rev  were  on  strong  diffusion 
(see  Fig  1.5,  Panel  (G),  22  Sep  1977,  0300-0500  UT)  indicating  longi¬ 
tudinal  extents  of  1  0000-1  5000  km. 

1 . 5  REP  -  Proton  Precipitation  Correlation 

In  scanning  through  the  data  set  it  was  immediately  apparent 
that  there  was  a  high  correlation  between  electron  and  proton  precipi¬ 
tation  (see  catalogue  remarks.  Appendix  C).  Particularly  noticeable 
was  the  similarity  between  precipitating  fluxes  in  the  235  and  435  keV 
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electron  and  the  80  keV  proton  channels.  An  example  of  this  feature  is 
presented  in  Fig  1.9  which  shows  the  12,  33,  235  keV  electron  channel 
and  80  keV  proton  channels  for  30  Jan  1977.  The  12  keV  electrons  are 
precipitating  from  19590  UT  to  20565  UT.  Precipitation  features  for 
33  and  235  keV  electrons  and  >80  keV  proton  are  evident  at  19646-19700 
UT,  19815  UT,  and  20475  UT  (proton  precipitation  began  at  19580  UT). 

This  suggests  the  precipitation  mechanism  for  both  electrons  and  protons 
may  be  either  related  or  identical.  Protons  did,  however,  precipitate 
much  more  frequently  than  electrons  and  their  zone  of  strong  diffusion 
precipitation  generally  extended  to  lower  latitude.  This  frequency  dif¬ 
ference  can  be  readily  seen  in  Fig  1.10  where  the  revs  (1-4  outer  zone 
passes  per  rev)  with  strong  electron  (Eprec.jp  >  235  keV)  and  proton 
(Eprecip  >  80  keV)  precipitation  normalized  to  the  number  of  revs  is 
plotted  for  increasingly  disturbed  magnetic  conditions.  Note  that  pro¬ 
ton  precipitation  is  seen  on  65%  of  revs  while  electrons  on  22%  of  revs 
for  all  magnetic  conditions.  The  high  occurrence  of  proton  precipita¬ 
tion  in  the  energy  range  (100-200  keV)  has  been  previously  observed 
(e.g.  On  ESR0  1A,  Linda!  en  et  al  . ,  1971). 

1 .6  Classification  of  REP  Events 

Due  to  the  strong  correlation  of  REPs  with  proton  (>80  keV) 
precipitation,  the  data  set  was  analyzed  for  REPs  without  concurrent 
strong  proton  precipitation.  Only  7  of  31 3  outer  zone  passes  fell  in 
this  category.  In  addition,  all  7  cases  had  concurrent  33  keV  elec¬ 
trons  precipitating  on  strong  diffusion  and  were  found  at  altitudes 
below  1 500  km. 
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FIGURE  1 . 


Flux  versus  (JT  plot  of  12,  33,  235  keV  electrons  and 
80  keV  protons  for  rev  1656  30  Jan  1977 
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FIGURE  1  .1 


0  Percentage  of  revs  (1-4  outer  zone  passes)  with  strong 
e  or  p  precipitation  normalized  to  revs  number  by 
magnetic  conditions. 
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Analysis  of  REP  energy  dependence  revealed  that  the  overwhelm¬ 
ing  majority  of  REPs  had  energies  in  the  range  235-1000  keV  on  strong 
diffusion.  Only  4  events  had  spectra  such  that  only  energies  above 
1  MeV  were  on  strong  diffusion.  These  same  4  cases  were  different  in 
another  way.  They  occurred  equatorward  of  the  outer  edge  of  the  outer 
zone;  the  remaining  cases  except  for  the  seven  above  all  occurred  right 
at  the  outer  edge  or  extended  inward  from  the  outer  edge. 

The  four  cases  with  strong  precipitation  only  above  1  MeV  were 
not  among  the  seven  cases  which  showed  no  concurrent  proton  precipita¬ 
tion.  REPs  could  thus  be  divided  into  three  distinct  classes:  those 
with  strong  diffusion  precipitation  below  (302)  and  only  above  (4)  1 
MeV  with  concurrent  protons  on  strong  diffusion  and  those  with  strong 
diffusion  precipitation  below  1  MeV  but  without  strong  proton  precipi¬ 
tation  (7).  These  classes  are  displayed  in  invariant  latitude  versus 
MLT  in  Fig  1.11  and  will  be  discussed  in  detail  in  Chapter  2. 


50 


FIGURE  1.11  REP  events  displayed  according  to  invariant  latitude 
and  MLT  for  the  following  classes:  dots  and  lines 
(.  I  )  for  strong  precipitation  at  energies  5  1  MeV 
with  concurrent  proton  preci pitation;  squares  (  o  ) 
for  strong  precipitation  at  energies  _>  1  MeV  with  con¬ 
current  proton  precipitation;  and  circles  (0)  for  strong 
precipitation  at  energies  <<  1  MeV  without  concurrent 
proton  precipitation. 
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CHAPTER  2.  THEORY  AND  CASE  STUDIES 


2.1  Background 

To  account  for  the  observed  subauroral  radio  signal  absorption, 

Bailey  (1968)  estimated  that  the  precipitating  electrons  at  energies 

3  2 

exceeding  several  hundred  keV  required  fluxes  Jp  >  10  electrons/cm  s 
str.  Such  intense  subauroral  (closed  field  lines)  precipitation  re¬ 
quires  pitch  angle  scattering  of  trapped  radiation  belt  electrons  near 
the  strong  diffusion  limit.  This  leads  to  an  isotropic  pitch  angle 
distribution,  i.e.,  precipitated  fl ux  trapped  flux  (Kennel ,  1969). 

The  rapid  scattering  can  occur  during  resonant  interaction  between 
plasma  waves  and  particles  in  the  magnetosphere.  Either  of  the  first 
two  adiabatic  invariants  (u  and  J)  can  be  violated  during  the  pitch 
angle  diffusion  and  this  can  provide  a  major  loss  mechanism  for  geo- 
magnetically  trapped  particles  (Thorne,  1975). 

The  condition  for  first  order  resonant  scattering  with  rela¬ 
tivistic  electrons  is: 

.  Q  .  SI. 

3,1  “  ujY 

where  =  "^r*  parallel  (to  ambient  magnetic  field) 

refractive  index 

Vo 

r»  -  c,  ,  normalized  parallel  electron  velocity 
<->  =  wave  frequency 

-*^=  relativistic  electron  gyrofrequency 

Y  =  (i  ~  j~>  )  ,  relativistic  mass 

enhancement  factor 
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c  =  speed  of  light  in  vacuo 
*ii  =  parallel  propagation  constant 
(wave  number) 

Thorne  (1977)  has  scaled  the  general  morphology  of  velocity  space  sur¬ 
faces  for  resonance  with  a  wave  of  a  given  frequency  and  refractive 
index  (Fig  2.1).  Waves  with  pn  >  1  generally  require  u>  <  Q  to  re¬ 
sonate  with  relativistic  electrons  because  only  this  surface  (^  <  1 ) 
reaches  the  loss  cone  at  relativistic  energies.  Waves  with  to  >  ft 
even  with  u,,  <  1  as  in  the  case  of  Auroral  Kilometric  Radiation  (Kurth 
et  a!  . ,  1975)  do  not  usually  resonate  with  relativistic  electrons.  A 
brief  review  of  plasma  waves  follows. 

Southwood  (1978)  has  reviewed  the  plasma  waves  found  in  the 
magnetosphere.  Four  distinct  classes  of  waves  have  been  identified 
in  the  outer  radiation  zone  4  <  L  <  8  where  REP  events  are  observed. 

The  first  class  are  electron  electromagnetic  waves  at  frequencies  be¬ 
low  the  electron  gyrofrequency;  waves  typically  occur  in  the  frequency 
band  .1  <  ~  <  .5.  These  so-called  right  hand  polarized  (Burton  and 
Holzer,  1974)  whistler  mode  "chorus"  waves  are  generated  at  high 
altitudes  (near  the  magnetic  equator)  outside  the  plasmasphere  with 
typical  wide  band  wave  amplitudes  JO  m  y  (Dunckel  and  Hell  iwell ,  1  969, 
Anderson  and  Maeda,  1977).  The  waves  propagate  down  the  field  lines  to 
low  altitudes  where  they  are  the  most  frequently  occurring  emission 
outside  the  plasmapause  (Kennel  and  Thorne,  1967;  Thorne  and  Kennel , 
1967;  Gurnett  et  al ■ ,  1969;  Dunckel  and  Hell iwell,  1969;  Tsurutani  and 
Smith,  1974).  Dunckel  and  Helliwell  (1969)  examined  emission 
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FIGURE  2.1  Velocity  space  surfaces  for  resonance  with  a  wave 
of  frequency  to  and  parallel  refractive  index  y  „ 
for  y  versus  y  „  . 
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intensities  for  300-500  Hz  whistler  mode  chorus  and  found  that  they 
peak  near  1000  MLT.  Also,  just  outside  the  plasmapause,  intensities 
increase  abruptly  at  0600  MLT  and  decrease  gradually  after  1400  MLT 
in  agreement  with  0V3-3  1.3  kHz  low  altitude  emissions  (McPherson  and 
Koons,  1970).  A  one  year  study  of  0G0-5  search  coil  magnetometer  data 
show  that  equatorial  chorus  is  strongly  substorm  correlated  with  two 
peaks:  one  at  post  midnight  (L=5-8)  and  a  second  in  the  dawn  to  noon 
sector  (L=7-ll )  (Tsurutani  and  Smith,  1977).  The  source  of  the  waves 
is  generation  by  loss  cone  instability  associated  with  the  pitch  angle 
anisotropy  of  substorm  injected  electrons  (Kennel  and  Petschek,  1966). 
Association  of  this  wave  phenomenon  with  precipitation  of  energetic 
electrons  (10-100  keV)  has  been  made  by  several  authors  ( Dungey,  1963; 
Cornwall ,  1964;  Kennel  and  Petschek,  1966;  Vampola  et  al ,  1971;  Thorne, 
1974;  Tsurutani  and  Smith,  1977). 

The  second  class  of  waves  are  electron  electrostatic  waves 
found  at  frequencies  above  the  electron  gyrofrequency  (Kennel  et  al . , 
1970).  These  so-called  "n  +  1/2"  waves  (^  z  n  +  1/2)  are  found 
throughout  the  magnetosphere  outside  the  plasmapause  (Kennel  et  al . , 
1970;  Shaw  and  Gurnett,  1975)  and  are  thought  to  be  generated  by 

>  0  electron  distributions  ( Fredricks ,  1971).  The  waves  are  lo¬ 
calized  about  the  geomagnetic  equator  from  L=4-10  a  majority  of  the 
time  in  the  morning  sector  ( Kennel  et  al . ,  1970;  Fredericks  and  Scarf, 
1973).  Wave  amplitudes  range  from  1-10  mV/m  ( Kennel  et  al . ,  1970) 
with  excursions  to  100  mV/m  during  substorms  (Scarf  et  al  . ,  1  973). 
Theoretical  parametric  studies  yielding  realistic  spatial  growth  rates 
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for  these  waves  have  been  accomplished  (Young  et  al  . ,  1  973;  Karpman 
et  al  . ,  1975;  Ashour-Abdalla  and  Kennel,  1976)  and  computations  of  the 
resonant  electron  population  show  energies  of  1-10  keV  (possibly  ac¬ 
counting  for  the  diffuse  aurora,  Lyons,  1974)  or  higher  (~  100  keV) 
during  substorms  (Scarf  et  al . ,  1973).  However,  this  class  of  waves 
does  not  resonate  effectively  with  relativistic  electrons  in  strong 
diffusion,  requiring  wave  amplitudes  in  excess  of  100  mV/m  to  cause 
strong  diffusion  of  100  keV  electrons  (see  Fig.  2.2 ,  Lyons,  1974). 

This  class  will  not  be  considered  further  in  this  study. 

Analagous  to  the  right  hand  electron  electromagnetic  waves 
which  cut  off  at  Q,  ,  there  is  a  third  class,  the  left  hand  ion  electro 
magnetic  waves  below  (e.g.,  Cornwall,  1965).  In  the  micropulsation 
literature  these  are  sometimes  referred  to  as  "PCI  or  IPDP  micropulsa¬ 
tions"  or  periodic  emissions:  PCI  ("Pearl s")are  dispersive  waves  pro¬ 
pagating  along  field  lines  thought  to  be  generated  by  the  loss  cone 
distribution  of  ring  current  protons  just  inside  the  plasmapause 
( Cornwal 1  ,  1  965;  Cornwall  et  al . ,  1970;  Mill  jams  and  Lyons,  1  974;  Kaye 
et  al . ,  1979);  "IPDP,"  intervals  of  pulsations  of  diminishing  period 
are  non-dispersive  waves  propagating  near  dusk  during  substorm  expan¬ 
sion  (McPherson  et  al . ,  1  968;  Lukkari  and  Kangas,  1  976;  Bossen  et  al . , 
1976).  In  a  statistical  study  by  Bossen  et  al .  (1976)  EMIC  occurred 

mostly  about  the  dusk  meridian  with  frequencies  ^  ~  .1-.2  and  wave 

+ 

amplitudes  l-7y  during  the  main  phase  of  geomagnetic  storms.  PCI  tur¬ 
bulence  has  also  been  observed  in  detached  plasma  regions  (Kivelson  et 
al . ,  1972)  prevalent  in  the  late  morning  to  dusk  sector  ( Chappel 1 , 
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GURE  2.2  Minimum  wave  amplitudes  required  for  electron  strong 

diffusion  scattering  at  L=7  versus  electron  energy  for 
electron  thermal  energies  0.1,  1  ,10  keV  (Lyons,  1  974  ) 
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Electron  Energy  (keV) 
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1974).  Theoretical  calculations  show  that  these  waves  are  candidates 
for  hard  REPs  (precipitating  electron  energies  >1  MeV)  in  the  evening 
plasmapause  (Thorne  and  Kernel,  1971)  and  possibly  during  the  day  in 
detached  plasmas  (Thorne,  1974). 

The  last  class  of  waves  to  be  discussed  here  are  the  electro¬ 
static  ion  cyclotron  waves  predicted  by  Kennel  et  al .  (1970)  after 
the  analogous  electrostatic  electron  wave  discovery  and  measured  by 
Gurnett  and  Frank  (1  977)  and  Kintner  et  al  . ,  (1  978).  These  waves 

have  frequencies  n+  <  tn  <  amplitudes  on  the  order  of  a  few 

/ 

10s  of  mV/m,  and  are  found  at  all  local  times  at  L  >  6  (Gurnett  and 
Frank,  1977;  Kintner  et  al . ,  1978).  Generation  of  the  waves  is  thought 
to  be  by  current  driven  instability  low  on  the  field  line  to  yield 
resonant  Landau  growth  (Kindel  and  Kennel,  1971;  Ashour-Abdall a  and 


Thorne,  1978)  and/or  equatorial  anisotropic  pitch  angle  instability 
due  to  distributions  of  hot  protons  (Coroniti  et  al . ,  1972;  Ashour- 
Abdall  a  and  Thorne,  1977;  Ashour-Abdall a  and  Thorne,  1978).  Observa¬ 
tions  have  indeed  correlated  REPs  with  intense  ELF  electrostatic  waves 
(Koons  et  al . ,  1972)  and  shown  that  some  cases  can  be  explained  only 
as  a  low  altitude  phenomenon  (Vampola ,  1977). 


2.2  Resonant  Populations 

Since  the  observational  data  covers  the  electron  energy  range 
12  keV-1 .6  MeV  and  proton  range  >  80  keV,  >  150  keV  etc.,  it  is  im¬ 
portant  to  determine  the  electron  and  proton  resonant  populations  for 
each  wave  phenomena  If  any  success  at  discrimination  (between  wave 
mechanism)  is  to  be  accomplished.  We  proceed  with  calculations  for 
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each  wave  type. 


2.2.1  Whistler  Mode  Chorus 

The  cold  plasma  dispersion  relation  for  parallel  whistler  mode 
chorus  propagation  assuming  «  w  «  Q  and  p  ,,  »  1  (cj^  >>  fi  )  is 
given  by  Kennel  and  Petschek  (1966): 

(2) 

where  w  =  plasma  frequency. 

Using  this  dispersion  relation  and  the  condition  for  first  order  cy¬ 
clotron  resonance  with  relativistic  electrons  (eq.  1)  assuming 
8„  -8  (valid  near  the  loss  cone),  Thorne  (1974)  has  derived  an 
expression  for  resonant  electron  energies: 
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,  magnetic  energy  per  particle 


8 tv  M* 

=  local  magnetic  field  value 

=  particle  number  density 
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=  nfc  ,  particle  rest  energy. 

w 

Thus  the  expression  for  the  resonant  electron  population  with  whistler 
mode  chorus  waves  of  frequency  oj  is: 
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where  E~  is  the  resonant  electron  energy. 
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B 

Since  B  :  -y  at  the  equator  where  Bo  =  .311  Gauss  and  N  varies  from 

L  3  3  3 

approximately  10  /cm  in  the  plasmasphere  ( L  =4 )  to  1-10/cm  in  the 

plasma  trough  region  ( L =5-8 )  for  a  K^  ~  2-3  (Chappel 1  et  al . ,  1  970), 

we  have  a  range  of  of  .5  keV  (1=4,  N  =1 0^ )  to  150  keV  (L=5,  N=1 )  at 

the  equator.  Down  on  the  field  line  (away  from  the  equator)  where  B 

increases  and  N  is  less  than  ionospheric  densities  {z  10^)  the  value 

4-3 

of  Em  can  reach  a  few  100  keV.  Values  of  N  not  greater  than  10  cm 
are  used  since  whistler  mode  chorus  propagating  down  the  field  line 
rapidly  becomes  perpendicular  to  the  local  magnetic  field  line  and 
reflects  at  frequencies  near  the  lower  hybrid  resonance  (Thorne  and 
Kennel ,  1  967).  [For  propagation  at  angles  to  the  field  line,  the 
cumulative  contribution  from  parallel  velocity  resonances  is  appropri¬ 
ate  since  higher  order  resonances  become  important  (Stix,  1962;  Kennel 
and  Petschek,  1966)  although  this  is  not  the  case  for  electromagnetic 
waves  with  kx  a  0.]  So  for  .1  <  .5  (Tsurutani  and  Smith,  1  977) 

and  =  .5-150  keV  (both  typical  equatorial  values),  we  have  E~  = 
1-800  keV.  Off  the  equator  (down  the  field  line  where  ^  is  lower 
and  Em  is  higher)  Er  is  expected  to  reach  values  of  nearly  10  keV. 

For  proton  precipitation  we  use  the  resonance  condition  for 
first  order  (n=l )  cyclotron  resonance  with  protons 


(5) 


and  the  dispersion  relation  for  whistlers  to  obtain 


63 


T  3  ^ 

where  —  =  1.837  x  10,  proton  to  electron  mass  ratio.  We  take  Em  ^ 
m" 

,  since  we  assume  the  densities  of  electrons  to  be  the  same  as  pro¬ 
tons  (DeForest  and  Mcllwain,  1971).  Using  the  same  range  of  values 
for  ~  we  can  readily  see  that  Ej!!  will  typically  be  greater  than  an 
MeV  (except  in  the  dense  pi asmasphere)  in  agreement  with  Kennel  and 
Petschek  (1 966) . 

Resonance  energies  for  both  electrons  and  protons  are  shown 
in  Fig.  2.3.  One  readily  sees  that  electrons  can  be  resonant  at 
energies  up  to  several  hundred  keV  while  protons  are  resonant  only 
at  energies  above  1  MeV. 


2.2.2  Electromagnetic  Ion  Cyclotron  ( EMI C )  Waves 

The  cold  plasma  dispersion  relation  for  parallel  propagating 
EMIC  waves  below  is  given  by  Kennel  and  Petschek  (1966): 


C) 


Thorne  and  Kennel  (1971)  and  Thorne  (1974)  have  used  equations  (1 ) 
and  (7)  to  derive  the  equation  for  resonant  electrons  with  EMIC: 
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FIGURE  2 . 


3  Electron  and  proton  resonant  energies  for  the  whistler 
mode  chorus  mechanism  (.1  <  £  <  .5)  [in  keV]  versus  L 
for  N=1 04cm~3  for  L=4  and  N=1  ,  10  cm“3  for  L>5. 
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and  we  see  that  (8)  becomes 


*•  3  2 

Using  the  previous  value  of  E~  at  L=4,  N=10  cm ",  and  ~  ^  .1-.5 
(Bossen  et  al . ,  1976)  we  see  that  E”  ~  few  MeV,  see  Fig  2.4.  Even 
during  intense  injection  events  where  EMI C  waves  may  be  generated 
outside  the  plasmasphere  to  L=8  due  to  enhanced  ion  fluxes  and  larger 
loss  cone  distributions  (Coroniti  et  al . .  1976),  E~  for  L=8,  N=10, 
and  ~  -  .2  is  5  MeV. 

The  equation  for  EMIC  resonance  with  protons  is  given  by 
Kennel  and  Petschek  ( 1 966 ) : 


Using  the  same  Em"  values  and  ^  =  .1-.2  (Bossen  et  al . ,  1976)  we 
see  E*  ~  1  Os-1 00s  keV  near  the  outer  edge  of  the  plasmasphere. 

In  summary,  EMIC  waves  resonate  with  electrons  only  at  ener¬ 
gies  >>  1  MeV  and  only  near  1  MeV  in  the  relatively  dense  plasmasphere 
or  detached  plasma  regions  where  N  is  relatively  large  and  at  the 
equator  where  B  is  lowest.  The  resonant  proton  population  for  the 
same  wave  mechanism  is  on  the  order  of  10's  and  100’s  keV. 

2.2.3  Electrostatic  Ion  Cyclotron  (ESIC)  Waves 

Since  the  frequency  of  ESIC  waves  is  near  the  I7+,  we  can 


FIGURE  2.4  Resonant  electron  energies  versus  L  for  the  ion 

cyclotron  {—  =  0.5)  wave  mechanism  for  densities 
+ 

_3+  3  _3 

N=1 ,  10  cm  outside  the  plasmasphere  and  N*1  0  cm 

L=4  (Thorne  and  Kennel ,  1971). 
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eliminate  the  unity  term  in  (1)  to  obtain: 


'S 


HD 


2  2  2 

Using  y  B  2.  y  -1  (valid  for  particles  near  the  loss  cone)  in  (11)  we 
have: 


*  *-  I 


Mx 


Si. 


(12) 


*  ,  2  „  2  r> 1 

2  C-  K»  -f.  St. 

By  writing  the  wave  refractive  index  yU„  =  ^j2  J  ,  where 

_ft-ris  the  ion  gyroratio,  is  the  ion  thermal  speed,  and  = 
Vth/n+.  ecl-  (12)  becomes: 


Si. 
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Vth,2  ^ 


Using  (~jrv  =  2  — where  E*h  is  the  ion  thermal  energy,  in  (13) 


we  have: 
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(14) 


The  resonant  electron 


energy  for  ESIC  waves  is  thus: 


(15) 


7C 


Ashour-Abdal 1  a  and  Thorne  (1977  and  1978)  have  solved  the  electrostatic 


dispersion  relation  numerically  and  obtained  k  space  contours  of  tem¬ 
poral  and  convective  growth  rates  normalized  to  +  for  typical  equa¬ 
torial  and  ionospheric  conditions.  If  waves  are  excited  by  the  current 
driven  instability  in  the  topside  ionosphere  with  electron  drifts  just 
above  the  Kindel  and  Kennel  (1971)  marginal  stability  limit,  1  Re 
above  the  earth  for  L  n  10,  E^  2.  2ev,  the  first  three  harmonics 
(typical  of  wave  observations,  Kintner  et  al . ,  1978)  have  typical 
values  for  kn  ^ r  —  .26-.1  (Ashour-Abdall a  and  Thorne,  1  978). 

This  yields  values  of  E^  -  50-300  keV.  For  excitation  in  the  equa¬ 
torial  plasmasheet  by  the  loss  cone  iuscabi  ity,  typical  equatorial 
parameters  are  E^  2  1  keV,  and  k„  cs  .4-. 5  for  the  first  five 
harmonics.  This  gives  2  1  00  keV  -  2  MeV.  Again  using  the  first 
order  cyclotron  relation  for  resonance  with  protons  (eq.  5)  we  have: 

si. - 

'  "/'v3-  '  ^ 


Rewriting,  we  have 


or 


Substituting  for 


we  have: 


•2 
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and 


_l_ 

L 


rrr^ 


vtt  rrr\  C. 

2  K.f  c* 


(17) 


thus 


(18) 


Since  (^—^— )2  (l-^)2  is  of  the  order  one  (1),  E*  -  E*h ,  i.e., 

2  ev  -  1  keV.  It  must  be  noted,  however,  that  higher  order  resonances 


must  be  taken  into  account  since  these  electrostatic  ion  cyclotron 


j. 

waves  have  polarizations  t—  -  1 0  (Kintner  et  al . ,  1978).  It  will  be 

K  ii 

shown  in  the  next  section  that,  unlike  electromagnetic  waves.,  resonant 
electrostatic  wave  amplitudes  required  for  strong  diffusion  scattering 
decrease  with  increase  in  energy  of  the  resonant  population  (y  -*•  ®) . 
Thus  resonant  energies  for  ions  can  easily  reach  a  few  100  keV 
(Ashour-Abdal la  and  Thorne,  1978).  So  in  summary,  electrons  with  E' 

-  50  keV-2MeV  and  ions  with  E*  -  1  »  100's  keV  can  resonate  with 

electrostatic  ion  cyclotron  waves. 


2 . 3  Wave  Ampl itudes 

Now  that  the  resonant  electron  and  proton  populatior 
known  for  each  wave  mechanism,  it  is  useful  to  determine  **• 
resonant  wave  amplitudes  required  for  strong  diffu'ic  n 
it  relate  the  ability  of  observed  wave  amplitude-,  * 


7^1  54b 

AIR  FORCE  INST  OF  TECH  *Rl SHT-PATTERSON 

AFB  OH 

relativistic  electron  precipitation:  an 

OBSERVATIONAL 

STUDY. (U) 

I960  L  J  ANDREOLI 

UNCLASSIFIED 

AFIT-CI-80-26D 

NL 

B 


■ 


H 


populations  on  strong  diffusion  but  the  equations  will  help  later  when 
examining  real  cases. 

The  condition  for  strong  diffusion,  i.e.,  very  rapid  scattering. 


it  ? 

C  'An 


where  =  diffusion  loss  time 
2^6=  1/4  bounce  time 


This  permits  the  particle  distribution  to  become  isotropic  even  over 

the  loss  cone  (=  >  J  ~  J+).  The  loss  time  approaches  a  minimum  value 

under  strong  diffusion  and  the  diffusion  coefficient,  D,  is  larger  than 

1/t.  ,  where  t  .  is  the  minimum  rate  of  diffusion  loss.  t.  ~ 
min  min  min 

-■^*3-  where  aL  is  the  pitch  angle  of  the  loss  cone. 

2.3.1  Electromagnetic  Waves 

The  action  of  electromagnetic  waves  on  particles  is  mainly  via 

pitch  angle  scattering.  The  pitch  angle  scattering  coefficient,  D  , 

cxcx 

is  given  by  Kennel  and  Petschek  (1966): 

*  (a*)1  -  it- Si-  M-V 

TTt  l  B-/  (20» 

where  b  =  wave  amplitude 

Bq  =  ambient  magnetic  field  value 
and  for  strong  diffusion  scattering  one  requires 


z  Tl 
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where  =  loss  cone  pitch  angle 

Dsd  =  coefficient  for  strong  diffusion. 

Equation  (20)  can  be  rewritten  including  the  relativistic  mass  correc¬ 
tion  and  bounce  averaging  to  obtain: 


where  f  =  is  the  fraction  of  the  electron  bounce  orbit  Spent  in 
resonance  with  waves. 

2  1 

For  equatorial  loss  cones,  a,  ~  — s-  .  For  electrons,  x,  MR  ; 
LR  LR  L  2l 

- -  -  where  R  =  1  earth  radius.  Equating  (22)  and  (20)  with 

v" 

substitutions  yields  a  fluctuating  field  wave  amplitude  for  strong 
electron  diffusion: 


3J  a  y  ^ 

Hi >C  R.  Jl.  f  ) 


(23) 


For  ions,  x,  /dR  =  )  n~  ,  n  =  —  fi+,  and  y~  1 

+  +  1/46  v+  (2E+/eV/2C  •  m  + 

where  E  and  EQ  are  the  Ion  energy  and  ion  rest  energy  respectively. 

Likewise,  the  fluctuating  field  wave  amplitude  for  strong  Ion  diffu¬ 
sion  is 
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and 
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Strong  electron  diffusion  (E~  ~  1/2  ■+  1  MeV;  y  2-3)  thus  re¬ 
quires  amplitudes  of  1  gamma  at  L=4  and  100  my  at  L=8.  For  ion  strong 
diffusion  (E^.~  50  keV),  amplitudes  of  3.5  gammas  at  L=4  and  300  my 
at  L=8  are  required.  These  values  are  within  the  range  of  amplitudes 
for  both  electrons  and  ions  for  EMIC  waves  (l-7y)  found  by  Bossen  et 
al .  (1976)  but  only  the  most  intense  chorus  events  with  amplitudes  ap¬ 
proximately  100  my  (Tsurutani  and  Smith,  1977)  meet  the  criteria  for 
electrons  at  L=3. 


2.3.2  Electrostatic  Amplitudes 

Electrostatic  waves  scatter  electrons  primarily  by  energy  dif¬ 
fusion.  Scattering  by  energy  In  the  perpendicular  (to  the  magnetic 
field  line)  direction  at  the  edge  of  the  loss  cone  requires  a  diffusion 
coefficient  as  follows: 
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For  strong  diffusion  which  can  maintain  isotropy  over  the  loss 
cone  one  requires 
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where  =  B (1 00  km). 

An  approximate  expression  for  Dj.*  as  a  function  of  electro¬ 
static  wave  amplitude  can  be  derived  (e.g.  Ashour-Abdalla  and  Thorne, 
1978): 
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where  p  is  the  electric  field  amplitude  of  the  electrostatic  wave. 
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Equating  eqs.  (27)  and  (28)  and  using  =  — —  ~  ,  we  have 

for  the  fluctuating  wave  amplitude  for  electron  strong  diffusion 
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For  y 


At  the  equator  pp  --W  and  £  =  [3.7  x  108  ^-t]1^2^). 

d.  o.  j  w«a  .  /  m 

£-  M  CU  m\f  *“ 

=20-250  V  for  L=8-4.  In  the  topside  iono- 
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sphere  near  1  Re,  f2-  a  To  ’  Bo  a  ,05,  and  t1  -4  x  10?  f]1/2^), 

For  *  =  2-3,  we  have  6.  ~  1-3  ^  for  L=8-4. 
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For  the  concomitant  Ion  strong  diffusion  using  eq.  (27)  and 
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we  have 
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Again  for  the  equator  _1_  andlT=l  we  have  £  =  6  - 

v  ^  D  i 

66  ^  at  L=8-4  for  a  50  keV  ring  current  proton.  At  the  1  Res  ^  ~  yjj 

B  ~  .05  and  yrz  1,  we  have  =  550  ~  -  750  ^p  for  a  50  keV  ring 

current  proton. 

Electrostatic  ion  cyclotron  waves  have  electric  field  ampli¬ 
tudes  typically  on  the  order  of  10-1 00' s  ^p  (Gurnett  and  Frank,  1977 
and  Kintner  et  al . ,  1978).  Required  wave  amplitudes  for  both  ions  and 
electrons  are  easily  met  by  these  waves  in  the  equatorial  regime  but 
not  in  the  topside  ionosphere. 

Table  2.1  lists  the  range  of  parameters  for  strong  ion  and 
electron  diffusion. 


2.4  Case  Studies 

An  example  of  each  of  the  classes  described  in  Section  1.6  is 
presented.  Each  example  was  selected  as  being  most  representative 
with  amount  of  peripheral  data  available  as  an  added  factor. 


2.4.1  Case  1,  ESIC  (rev  479,  day  250,  6  SEP  1976). 
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TABLE  2.1  Theoretical  REP  mechanisms  and  characteristics. 
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and  protons  80-770  keV  for  three  outer  zone  passes  are  shown  in  Fig  2.5 

а, b,c,  respectively.  The  low  altitude  outer  zone  pass  at  local  evening 
(Fig  2.5a)  shows  an  extended  region  of  12  keV  precipitation  (isotropic) 
at  L  »  5-10.  This  is  actually  an  inverted  V  event  which  is  bounded  by 
235  keV  precipitation.  At  L=5  the  precipitation  extends  to  nearly  1  MeV 
electrons  and  80  and  150  keV  protons.  This  agrees  well  with  the  cri¬ 
teria  of  ESIC  (Table  2.1).  Fig  2.5b  shows  the  outer  zone  pass  at  high 
altitude  (~  8000  km)  near  local  dusk.  Twelve  (12)  keV  electrons  are 
precipitating  from  L  *  7-9  with  two  large  spikes  near  L  ~  8.4.  The 

two  spikes  are  associated  with  two  electrostatic  shocks  (Kintner,  pri¬ 
vate  communication)  with  perpendicular  electric  field  £x  *  500  mV/m. 

In  addition  >  80  and  >  150  keV  protons  are  precipitating  from  L  £  6. 
Although  there  is  significant  trapped  235  keV  electron  flux  from  L 

б.  8-7. 8,  there  is  no  precipitation.  The  next  pass  (Fig  2.5c)  is  a 
medium  altitude  (~  5000  km)  pass  near  local  morning.  Twelve  (12)  keV 
electron  precipitation  is  evident  from  L  a  6-13.  Extensive  33  keV  is 
coincident  although  not  on  strong  diffusion  from  L  ~  9-12.  The  235 
keV  channel  shows  "moderate"  precipitation  at  L  »  6-8.  Significant 
fluxes  of  >  80  keV  protons  are  present  to  roughly  L  *  12  without  any 
precipitation.  This  fits  the  criteria  of  whistler  mode  chorus  in 
Table  2.1. 

Pitch  angle  plots  of  the  low  altitude  local  evening  outer  zone 
pass  of  rev  479  are  shown  In  Fig  2.6.  Panel  (A)  at  5962  UT  shows  a 
normal  pitch  angle  distribution  both  down  (0°)  and  up  (180°)  the  local 
magnetic  field  line.  Panel  (B)  shows  the  next  distribution  towards 
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FIGURE  2.5a  Flux  (counts/sec,  every  fourth  data  point)  versus 

L  value  of  low  altitude  (700  km)  outer  zone  pass  near 
local  evening.  Vertical  arrows  (+)  indicate  singular 
precipitation  and  horizontal  arrows  (-*  <-)  indicate 
limits  of  extended  precipitation. 


81 


FIGURE  2.5b  Same  plots  for  high  altitude  (8000  km) 
local  dusk. 
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FIGURE  2.6  Pitch  angle  versus  flux  (counts/sec)  for  235-1610  keV 
electron  channels  for  3  times  in  the  local  evening 
outer  zone  of  rev  479.  The  upward  loss  cone  (L,C.  =  50° 
is  indicated  by  the  vertical  line. 


the  outer  edge  of  the  outer  zone  with  a  large  spike  In  the  upward  view¬ 
ing  loss  cone.  The  very  next  distribution  Panel  (C)  shows  Isotropic 
flux  in  the  upward  loss  cone.  Isotropic  fluxes  are  evident  to  ~  1  MeV. 

Since  the  precipitating  flux  Is  higher  than  the  locally  trapped 
flux,  one  can  estimate  the  altitude  at  which  pitch  angle  scattering 
occurred  (Vampol a,  1977).  This  is  accomplished  by  tracing  the  preci¬ 
pitation  flux  to  the  shape  of  the  outer  zone  as  illustrated  In  Fig  2.7. 

4 

At  L=4  the  trapped  flux  would  be  roughly  2x10  counts/sec  up  the 
field  line  where  the  precipitation  occurred.  Fig  2.5b  shows  fluxes 

4 

near  L=4  of  approximately  2x10  counts/sec.  The  altitude  of  this 
pass  is  ~  7600  km.  It  must  be  noted  that  this  "spikey"  event  is  very 
rare  with  the  overwhelming  majority  of  precipitation  cases  showing  flat 
isotropic  fluxes. 

Resonance  with  1  MeV  electrons  at  ~1  Re  would  require  -£ 

~  .008,  slightly  smaller  than  values  for  first  three  harmonics  calcu¬ 
lated  by  Ashour-Abdalla  and  Thorne  (1978).  In  addition,  the  wave  ampl  1 

tude  is  <f  *  350  — .  Both  these  calculations  Indicate  that  higher  har 
w  m 

monies  may  play  a  very  Important  role  In  the  mechanism. 

Fig  2.8a  shows  time  spectrograms  for  the  filtered  plasma  wave 
data  (frequency  .03-100  KHz)  and  low  energy  electron  (.2-33  keV)  and 
protons  (.1-3.9  keV),  235  keV  electrons  and  >  80  keV  protons  of  the 
low  altitude  local  evening  pass.  The  80  keV  proton  and  235  keV  elec¬ 
trons  are  seen  to  precipitate  (upward  loss  cone  filled)  at  5815  UT  and 
5950  UT  (the  spike  occurred  here).  Between  these  times  Is  the  Inverted 
V  electron  event.  Plasma  turbulence  Is  also  evident  but  not  well 
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FIGURE  2.7  Enlarg«aent  of  portions  of  Figure  2.5a. 
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FIGURE  2.8a  Spectrogram  of  plasma  waves,  235  keV  electrons,  low 
energy  electrons,  low  energy  protons,  80  keV  protons 
for  the  low  altitude  local  evening  pass. 
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defined.  Fig  2.8b  shows  the  high  altitude  local  dusk  pass  with  the  in¬ 
verted  V  electrons  evident  from  10500  UT  to  10700  UT  with  ion  beams 
also  evident  between  these  times.  Plasma  turbulence  is  also  evident 
between  these  times  but  the  235  keV  electrons  have  been  depleted. 

Fig  2.9  shows  the  raw  magnetogram  for  Kiruna  (67°  50'N,  20°  25'E) 
for  4-6  September  1976.  The  period  is  relatively  quiet  from  0200  UT  to 
2200  UT  5  September  1976.  A  small  substorm  occurs  shortly  before  mid¬ 
night  on  6  September  and  the  estimate  (to  be  shown  later)  of  the  extent 
of  the  precipitation  in  time  is  indicated  by  the  arrows.  In  addition, 
the  geosynchronous  satellite,  ATS-6,  showed  two  small  excursions  in  the 
140-800  keV  electron  channel  (see  Fig  2.10). 

The  S3-2  satellite  has  a  magnetic  spectrometer  similar  to  S3-3 
(but  only  to  electron  energies  of  400  keV).  These  satellites  were  co- 
planar  at  this  time  of  year.  The  satellite  tracks  are  shown  in  Fig  2.11. 
The  X's  indicate  regions  of  precipitation  for  either  satellite.  The 
precipitation  is  first  seen  on  S3-3  at  1.6  UT  at  21.5  MLT  (spike  event). 

A  similar  spike  event  is  then  seen  at  2.1  UT  (both  200  and  400  keV 
channels)  by  S3-2  in  the  same  location  and  altitude.  Neither  satellite 
observed  precipitation  electrons  (protons  were  precipitating.  Fig  2.5b) 
at  local  dusk.  Both  satellites  observed  precipitation  near  200  keV  at 
local  morning  at  3.05  UT  (S3-2)  and  3.8  UT  (S3-3).  The  S3-2  satellite 
continues  in  orbit  till  precipitation  at  200  keV  is  observed  again  at 
3.65  UT  at  3.0  MLT,  3.7  UT  at  2.3  MLT,  and  3.8  UT  at  22.5  MLT.  The 
extent  of  precipitation  is  thus  c2.2  hrs.  The  actual  S3-2  flux  versus 
times  are  shown  in  Fig  2.12a,b,c. 
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FIGURE  2.9  Kiruna  raw  magnetogram  for  4-6  September  1976. 
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FIGURE  2.10  ATS-6  for  6  September  1976.  Arrows  indicate  beginning 
and  ending  of  event. 
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FIGURE  2.11  S3-2  (1,  2,  3)  and  S3-3  satellite  (a,  b,  c)  tracks  for 

6  September  1976  in  invariant  latitude  versus  MLT. 
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FIGURE  2.12a  S3-2  flux  versus  time  plots  of  pass  1,  Figure  2.11. 


Arrows  Indicate  precipitation. 
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FIGURE  2.12c  Same  for  pass  3,  Figure  2.11. 
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In  summary  this  event  fits  the  criteria  for  ESIC  mechanism  at 
local  evening  with  whistler  mode  chorus  at  local  morning.  This  par¬ 
ticular  event  seems  to  be  a  low  altitude  (~  1.2  Rq)  phenomenon  corre¬ 
lated  with  a  1 20y  substorm  and  inverted  V  event. 

2.4.2  Case  2,  EMIC  (rev  2484,  day  132,  12  May  1977) 

The  flux  versus  time  plot  for  the  EMIC  example  is  shown  in 
Fig  2.13a,b.  No  precipitation  is  seen  in  the  12  or  33  keV  electron 
channels  for  this  local  evening,  high  altitude  ( ~  7800  km)  pass.  In 
the  235  keV  channel,  however,  the  flux  at  L  ~  4.3  falls  an  order  of 
magnitude  into  the  loss  cone  where  normal  loss  cones  are  nearly  2 
orders  of  magnitude  deep.  Proceeding  up  in  energy,  the  loss  cone  at 
this  same  L  value  begins  to  fill  until  at  850  keV,  the  loss  cones  are 
nearly  full  (Fig  2.13b).  Isotropic  distributions  are  evident  in  the 
1.1,  1.3,  and  1.6  MeV  channels.  In  addition,  the  >  80  keV  protons  are 
on  strong  diffusion  from  just  equatorward  of  the  electron  precipitation 
and  continuing  poleward.  The  loss  cone  filling  with  increasing  energy 
is  evident  in  Fig  2.14  where  pitch  angle  plots  of  flux  for  160  keV  to 
1600  keV  electron  energies  are  shown.  The  upward  observing  loss  cone 
is  at  0°  in  these  plots.  This  spectrum  and  concurrent  proton  preci¬ 
pitation  is  indicative  of  the  EMIC  mechanism. 

The  spectrogram  for  this  pass  is  shown  in  Fig  2.15.  Low  energy 
electron  and  proton  activity  is  evident  poleward  of  the  REP  and  only 
the  >  80  keV  is  coincident.  No  plasma  data  was  available  for  this 
pass. 

No  direct  means  of  observing  the  plasmapause  was  available  for 
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FIGURE  2.13b  Same  for  higher  energy  electrons  and  protons. 
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figure 


14  Pitch  angle  versus  flux  plots  for  electrons  with 
energies  160-1600  keV  for  12  May  1977. 
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this  pass.  However  Mizera  (1974)  has  shown  that  the  midnight  plasma- 


pause  is  ~.8  L  equatorward  of  the  energetic  proton  precipitation  boun¬ 
dary.  The  electron  and  proton  precipitation  are  nearly  coincident  at 
local  evening.  Since  the  plasmasphere  bulges  poleward  from  midnight 
to  dusk,  the  REP  may  indeed  be  just  within  the  plasmapause.  Another 
rough  approximation  to  the  L  value  of  the  midnight-dawn  plasmapause 
location  is  given  empirically  by  Carpenter  and  Park  (1973)  from  ATS-5 
data  : 


Lpp  =  5.7  -  0.47  Kp  (29) 

where  Kp  is  the  maximum  geomagnetic  planetary  three-hour-range  indices 
12  hours  previous  to  the  event.  The  same  formula  may  be  used  with  Kp 
as  the  actual  value  at  the  time  of  the  event.  For  UT  of  4.6  hrs  on 
12  May  1977  the  Kp  at  the  time  of  event  was  ~  4.5  and  the  max  12  hrs 
previous  was  5.0.  These  two  values  yield  Lpp  =  3. 3-3. 9  as  the  mid¬ 
night-dawn-noon  plasmapause  and  approximately  twice  that  value  for 
the  dusk  plasmapause.  These  are  sketched  in  Fig  2.16.  For  these 
rough  approximations,  the  location  of  the  precipitation  of  rev  2484 
is  clearly  within  the  plasmasphere. 

Finally,  Fig  2.17  shows  the  ATS-6  electron  data  on  12  May 
1977.  The  REP  occurred  during  a  highly  disturbed  period  as  indicated 
in  the  140-600  keV  electron  channel.  In  addition,  S3-3  showed  preci¬ 
pitation  (ESIC)  at  energies  up  to  435  keV  on  the  preceding  and  follow¬ 
ing  passes  at  the  same  local  time.  Only  four  of  these  events  (Ep  _> 

1  MeV)  were  observed  for  the  entire  data  set. 
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FIGURE  2.16  Plasmapause  relation  to  rev  2484, 
on  rough  estimates. 
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2.4.3  Case  3,  Chorus  (rev  1413,  day  366,  31  December  1976) 

The  flux  versus  time  plot  for  31  December  1  976  is  shown  in 
Fig  2.18.  Only  33  and  235  keV  electron  precipitation  is  evident  just 
prior  to  local  noon  at  low  altitude  (550  km).  No  protons  were  preci¬ 
pitating  on  strong  diffusion  though  they  were  present.  The  pitch 
angle  plot  for  this  event  is  shown  in  Fig  2.19  where  the  loss  cone 
55°.  Normal  distributions  are  evident  at  12812  UT  but  at  the  later 
times,  33  and  235  keV  channels  are  isotropic  in  the  upward  viewing 
loss  cone. 

The  spectrogram  (Fig  2.20)  shows  evidence  of  plasma  noise 
just  above  30  Hz  at  and  poleward  of  the  REP  (12820  UT).  A  faint 
band  near  10  kHz  is  also  evident. 
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FIGURE  2.1 


9  Pitch  angle  distributions  for  33,  235,  435  keV  electron 
channels  for  31  December  1  976. 
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FIGURE  2.20  Spectrogram  for  low  altitude  early  morning  pass 
on  31  December  1  976. 
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CHAPTER  3.  ATMOSPHERIC  CONSEQUENCES 


3.1  Introduction  and  Background 

The  earth's  atmosphere  below  100  km  has  been  the  object  of 
intense  investigation  over  the  past  10-15  years.  It  is  well  known 
that  a  variety  of  electromagnetic  and  particle  radiation  can  deposit 
energy  in  this  region  and  thus  modify  the  local  chemistry.  The  re¬ 
mainder  of  this  section  describes  the  ionization  and  neutral  chemistry 
with  special  emphasis  on  ozone  removal  .  In  subsequent  sections  we 
examine  the  three  examples  in  Chapter  2  and  their  effect  on  ozone. 

The  ionosphere  below  100  km,  known  as  the  "D  region,"  has  been 
studied  by  several  authors  with  respect  to  ionization  by  solar  and 
galactic  cosmic  rays  (Webber,  1962),  primarily  in  the  polar  cap  region 
(A  >  75°) :  energetic  electron  precipitation  in  the  auroral  region, 

Al  70-75°  (Rees,  1963;  Lerfold  and  Little,  1964;  Bail ey  et  al  . , 

1970)  and  subauroral  regions,  A  z.  60-70°  (Bail ey,  1968;  Larsen,  1973; 
Reagan,  1977;  Thorne,  1977b)  and  at  mid  latitudes,  A  <  60°  (Potemra 
and  Zmuda,  1970;  Spjeldvik  and  Thorne,  1975a,b;  Larsen  et  al . ,  1976). 
Solar  UV,  especially  HLya  at  1216  A°  and  HLyB  at  1025  A°,  produces  the 
dominant  primary  ion  in  the  quiet  mesosphere  through  ionization  of  NO. 
Also,  solar  UV  scattered  from  the  sunlit  hemisphere  has  long  been 
thought  to  be  an  important  ionizing  radiation  source  at  night.  Solar 
UV  (and  X-ray)  radiation  is  strongly  absorbed  below  about  85  km  where 
positive  ion  composition  is  dominated  by  multiply  hydrated  ion  clusters 
of  the  type  H+(H20)n,  N0+(H20)n,  and  Hj0+(H20)n  (Johannesen  and 
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Krankow  ski,  1972).  0^,  the  dominant  negative  ion  below  80  km,  is 

produced  by  electron  attachment  and  interacts  with  CO,  NO,  HCO  and 
(H20)n  to  form  heavy  negative  ion  and  large  water  clusters. 

Although  the  precipitating  particle  energy  is  usually  negli¬ 
gible  compared  to  other  ionizing  sources  below  100  km,  it  can  be  the 
dominant  ionization  source  during  disturbed  conditions  below  80  km 
(Thorne,  1977b).  Particulate  energy  produces  mostly  primary  ions 
0*  and  N^  ,  N+,  and  0+  of  which  -94%  are  rapidly  converted  to  0^ 
through  charge  exchange  with  O2  (Reagan  et  a!.,  1978)  and  -6%  to  N0+ 
by  ion-atom  interchange  with  ^  and  O2  (Gunton  et  al  . ,  1  977).  This 
disturbed  increase  in  0^  can  thus  modify  the  complex  clustering  reac¬ 
tions  in  the  D-region. 

A  comparison  of  several  ionizing  sources  is  presented  in  Fig 
A. 15  (Appendix  A)  where  the  overwhelming  dominance  of  the  largest 
solar  proton  event  (SPE)  of  August  1972  is  evident.  Immense  SPEs  such 
as  August  1972  are  relatively  rare  and  indeed  REPs  can  frequently  be 
the  dominant  ionization  source  below  100  km  (Spjeldvik  and  Thorne, 

1  975a;  Thorne,  1  977b). 

The  production  and  removal  of  odd  oxygen  (0,  03)  is  of  para¬ 
mount  importance  in  neutral  atmospheric  chemistry  below  100  km.  Odd 
oxygen  is  produced  primarily  by  photodissociation  of  molecular  oxygen. 
The  Chapman  (1  930)  process,  0+0^  -*■  2O2,  is  now  well  known  to  account 
for  only  -20%  removal  of  odd  oxygen.  The  major  removal  process  in¬ 
cludes  catalytic  loss  cycles  with  chlorine  compounds  (e.g.  Crutzen 
and  Howard,  1978),  and  odd  nitrogen  and  odd  hydrogen  species  (Blake 
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and  Lindzen,  1973).  The  latter  two  are  of  particular  interest  here 


since  they  can  be  significantly  modified  by  intense  particle  precipi¬ 
tation  (Thorne,  1977a, b;  1979). 

The  bulk  of  odd  oxygen  depletion  is  through  the  following 
nitrogen  chemistry  important  particularly  below  50  km  (Blake  and 
Lindzen,  1973): 


NO  +  03  -*  N02  +  02  (Davis  et  al  , ,  1973) 

N02  +  0  ■+  NO  +  02  (Herron  and  Huie,  1  972) 

REPs  can  be  the  dominant  source  of  NO  production  by  dissociation  of  N2 
in  the  mesosphere  and  a  competing  source  with  nitrous  oxide  upwelling 
from  the  Earth's  surface  in  the  upper  stratosphere  (Thorne,  1979). 

Above  40  km  the  odd  hydrogen  chemistry  begins  to  dominate 
In  the  removal  of  odd  oxygen  (Blake  and  Lindzen,  1973).  Here  the  odd 
hydrogen  is  produced  mainly  by  dissociation  of  water  vapor  by  photons 
and  0  ( ' D)  atoms  and  as  a  by-product  of  02  production  by  energetic 
particle  precipitation  (Rowe  et  al . ,  1974).  In  this  way,  the  increase 
in  02  during  REPs  may  also  affect  the  mesospheric  removal  of  odd  oxy¬ 
gen.  Since  energy  deposition  by  REPs  can  dominate  in  the  mesosphere, 
we  shall  subsequently  examine  the  removal  of  odd  oxygen  by  the  en¬ 
hanced  odd  hydrogen  chemistry  and  its  ozone  consequences  for  the  three 
REP  examples  of  Chapter  2. 

3.2  REP  Energy  Spectra  and  Atmospheric  Energy  Deposition 

The  precipitation  spectra  for  the  three  cases  presented  in 
Chapter  2  are  shown  in  Fig  3.1.  The  ESIC  spectrum  (Case  1)  has  the 
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FIGURE  3.1  Energy  spectra  for  cases  1,2,3  of  Chapter  2. 
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highest  flux  in  the  100-800  keV  range.  The  EMIC  spectrum  (Case  2) 
dominates  at  energies  above  1  MeV  and  is  comparable  to  the  previously 
published  spectrum  of  Vampol a  (1971).  The  chorus  event  is  relatively 
weak  dominating  only  near  3  and  30  keV.  These  REP  spectra  were  input 
to  the  UCLA  (Spjeldvik,  1974)  electron  precipitation  model  (see  Ap¬ 
pendix  D).  The  ion  pair  production  (assuming  35  keV  per  pair)  versus 
altitude  is  shown  in  Fig  3.2.  The  Bremsstrahl ung  X-ray  deposition 
profiles  were  also  plotted  in  Fig  3.2.  As  expected,  the  ESIC  curve 
dominates  in  the  70  km  region,  and  EMIC  in  the  50-60  km  region. 

Whistler  events  are  relatively  non-competitive  with  other  REP  spectra 
but  may  dominate  local  chemistry. 

The  ESIC  case  of  Chapter  2  was  used  to  model  a  typical  REP  to 
calculate  the  total  energy  deposition.  The  latitudinal  extent  was 
approximately  200-300  km  based  on  the  invariant  latitude  excursions 
of  both  S3-2  and  S3-3.  The  longitudinal  extent  based  on  local  time 
duration  was  7,000-10,000  km.  The  event  lasted  approximately  2.2  hours. 

The  peak  ion  production  rate  occured  at  ~73  km  altitude  with  a  value 
-  -3  -1 

of  6.8  x  10  (cm"  sec  ).  Assuming  precipitation  was  nearly  constant 

and  uniform  throughout  the  events  (which  is  not  likely),  a  peak  value 

24  -1  22 

of  ~1  x  10  ion  pairs  cm  (column)  or  3.5  x  10  keV  of  energy  was 

deposited  during  this  event. 

3.3  REP  Induced  Hydrogen  and  Ozone  Modifications 

The  REP  ion  production  rates  versus  altitude  were  processed 
to  produce  odd  hydrogen  (and  thus  hydroxyl)  according  to  the  chemical 
modelling  of  Heaps  (1978).  Fig  3.3  presents  this  model  of  odd  hydrogen 
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FIGURE  3.2  Ion  pair  production  rate  versus  altitude  for  ESIC, 
EMIC,  and  chorus  events. 
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FIGURE  3. 


Altitude  profile  of  the  odd  hydrogen  formation  rate 
related  to  ion  pair  production. 
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production  as  a  function  of  ion  pair  production  for  disturbed  condi¬ 
tions  (higher  ion  production  rate).  The  higher  altitude  limit  for 
cluster  ion  formation  is  clearly  indicated  by  the  trailoff  above  80  km. 
Since  the  nitrogen  cycle  dominates  in  the  stratosphere,  odd  hydrogen 
is  not  expected  to  be  important  below  40  km. 

After  modelling  the  OH  increases  to  the  REP  ionization  rates, 
a  numerical  simulation  fashioned  after  Blake  and  Lindzen  (1973)  was 
performed  to  calculate  the  relative  change  in  ozone  concentrations. 

The  numerical  simulation  differs  from  Blake  and  Lindzen's  model  in  that 
two  reactions  (OH  nighttime  source  reaction  and  Chapman's  (1930)  orig¬ 
inal  3  body  atomic  oxygen  collision  equation  to  form  0^)  are  added, 
reaction  rates  are  updated,  and  solar  flux  and  o zone  conentrations 
vary  with  time  of  day  (DeVore,  1977). 

The  results  for  Case  1  and  2  (Fig  3.4a  and  b)  are  similar  since 
they  both  occur  during  local  evening.  Factors  of  10  and  4  increase  in 
OH  concentr  *:ions  at  70  km  immediately  after  precipitation  commenced 
occurred  for  the  ESIC  and  EMIC  cases  respectively.  The  enhancement 
continued  throughout  the  night  and  several  hours  into  the  day.  Ozone 
was  not  affected  until  sunrise  due  to  the  low  nighttime  0  concentra¬ 
tion  but  showed  sunrise  decreases  of  10-15%  for  both  cases.  These 
diminished  ozone  concentrations  lasted  for  several  hours. 

The  chorus  event  (Case  3,  Fig.  3.4c)  occurred  near  local  noon 
with  a  small  enhancement  in  OH  concentration  and  throughout  the  suc¬ 
ceeding  night.  The  ozone  concentration  showed  an  immediate  response 
dropping  by  nearly  20%.  The  diminution  of  ozone  ended  at  local  sunset. 
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FIGURE  3.4a,b,c  Ozone  and  hydroxyl  concentrations  [cm  ]  versus 
local  time  (hrs)  for  a.  ESIC,  b.  EMIC,  and  c. 
chorus  REP  events. 
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UNPERTURBED 


A  similar  simulation  with  an  intense  REP  modelled  after  the 
ESIC  case,  this  time  occurring  at  local  noon,  showed  only  a  17%  in¬ 
crease  in  OH  but  more  than  a  30%  decrease  in  ozone. 

Altitude  profiles  of  the  ozone  concentrations  at  the  local 
time  of  maximum  decrease,  OH  concentrations  at  the  same  time  and  at 
time  of  maximum  OH  enhancement  for  the  three  cases  are  presented  in 
Figs  3.5a,b,c.  Maximum  ozone  decreases  occurred  at  69  km  for  both 
the  ESIC  and  EMIC  cases  and  at  76  km  for  the  chorus  case.  Only  in  the 
chorus  case,  where  onset  of  OH  enhancement  and  ozone  diminution  were 
simulatneous,  did  the  deviations  occur  at  the  same  altitude. 
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FIGURE  3.5a 


b,c  Altitude  [km]  profiles  of  ozone  concentration  at 
time  of  maximum  decrease,  hydroxyl  concentration 
at  the  same  time,  and  hydroxyl  concentration  at 
the  time  of  maximum  enhancement  for  a.  ESIC,  b. 
EMIC,  c.  chorus  REP  events. 
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CHAPTER  4.  SUMMARY,  CONCLUSIONS  AND  OUTLOOK 


I  have  presented  here  the  results  of  two  years  of  data  analysis 
of  fourteen  months  of  satellite  outer  radiation  zone  coverage  (Chapter 
1).  Theory  and  examples  have  been  provided  in  an  attempt  to  explain 
the  types  of  REPs  discovered  in  the  analysis  (Chapter  2).  Finally, 
some  atmospheric  consequences  were  investigated  (Chapter  3). 

For  the  first  time  the  properties  of  strong  diffusion  REPs  have 
been  studied  using  a  relatively  extensive,  dense  satellite  data  set 
and  the  following  results  have  been  obtained: 

o  REPs  show  a  seasonal  variation  with  equinoxial  maxima 
and  solstice  minima  in  agreement  with  Bailey  (1968). 
o  The  occurrence  frequency  of  REPs  is  nearly  linear  with 

general  magnetic  conditions  again  in  agreement  with  Bail ey 
(1968). 

o  Correlation  of  REP  producing  substorms  and  substorm  in¬ 
dices  (AE)  are  weak  in  the  limited  number  of  cases  exam¬ 
ined  in  agreement  with  Thorne  and  Larson  (1976)  and  not 
all  substorms  produce  REPs  in  agreement  with  Rosenberg 
et  al .  (1972)  and  Thorne  and  Larson  (1976). 
o  Overall,  only  one  out  of  ten  REPs  occur  in  daytime  but 
as  magnetic  conditions  become  disturbed  there  is  an  in¬ 
crease  in  daytime  REPs  approaching  three  out  of  ten. 
o  REPs  are  estimated  to  extend  in  longitude  from  1000  km 
to  10-15,000  km  in  general  agreement  with  Larsen,  1973, 
and  Bailey,  1 966. 
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o  98 %  of  all  REPs  occur  concurrently  with  protons  [2.  80  keV) 

precipitating  on  strong  diffusion  and  many  features  of  the 
precipitation  are  similar. 

0  Three  distinct  classes  of  REPs  occurred:  an  overwhelming 
majority  (96.5%)  had  strong  diffusion  precipitating  spectra 
from  hundreds  to  a  thousand  keV  energy  and  all  these  events 
occurred  near  the  outer  bound  of  trapping  simultaneously 
with  proton  precipitation;  the  smallest  group  (1.3%)  had 
strong  diffusion  precipitating  spectra  only  for  electrons 
>1  MeV  and  again  with  concurrent  proton  precipitation  but 
inward  from  the  outer  bound  of  trapping  and  probably  within 
the  pi asmasphere;  a  slightly  larger  group  (2.2%)  had  strong 
diffusion  precipitation  spectra  only  for  electrons  of  a 
few  hundred  keV  but  with  protons  present  with  normal  loss 
cone  distributions  (not  precipitating  on  strong  diffusion). 

The  theoretical  investigation  shows  that  REPs  can  be  caused  by 
radiation  belt  plasma  interactions  with  plasma  waves  and  that  the  elec¬ 
trostatic  ion  cyclotron  (ESIC),  electromagnetic  ion  cyclotron  ( EMI C ) 
and  whistler  mode  chorus  waves  can  respectively  account  for  the  three 
distinct  types  of  REPs  discussed  above. 

Atmospheric  consequences  were  investigated  with  the  following 

resul ts : 

0  ESIC  and  EMIC  REPs  deposit  more  energy  into  the  atmosphere 
than  chorus  REPs. 

o  Immediate  increases  in  the  OH  concentration  up  to  a  factor 
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of  10  are  observed  during  REP  occurrence  and  gradually 
return  to  normal  within  24  hours. 

o  Ozone  diminution  only  occurs  simultaneously  with  daytime 
REPs,  only  after  sunrise  with  nighttime  REPs,  and  return 
to  normal  in  several  hours. 

The  following  main  conclusions  may  be  drawn  from  these  results: 

o  The  ESIC  wave  interaction  is  by  far  the  dominant  cause  of 
REPs. 

o  REPs  are  locally  important  in  reducing  ozone  concentration 
(10-30%). 

Based  on  all  of  the  above,  future  considerations  should  in¬ 
clude  the  following  areas  of  endeavor: 

o  Although  the  concurrent  proton  precipitation  seems  to 
point  to  the  anisotropic  pitch  angle  instability  at  the 
equator  rather  than  to  a  current  driven  instability  low 
on  the  field  line,  results  are  inconclusive  (Kintner  et  al , 
1979).  Resolution  of  the  location  of  the  instability 
would  lead  to  better  understanding  of  the  trigger  mechanism 
for  REPs. 

o  Relativistic  electron  precipitation  at  rates  less  than 
strong  diffusion  were  not  analyzed  here  due  to  the  diffi¬ 
culty  with  distinction  of  the  pitch  angle  distribution  and 
the  Instrument  response  limits  for  small  loss  cones. 

Weaker  events  (moderate  diffusion)  may  be  locally  im¬ 
portant  to  ozone  chemistry  particularly  if  they  occur 
during  the  daytime. 
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o  Though  REPs  locally  increase  the  OH  concentrations  day  or 
night,  an  investigation  of  long  term  effects  on  the  atmo¬ 
sphere  are  important  since  it  may  be  possible  to  dry  out 
the  mesosphere  over  a  year's  time  and  actually  account  for 
long  term  ozone  buildups  ( Crutzen,  private  communication, 
1980). 
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Appendix  A  -  OV 1-19  Data  Set  and  Analysis 

The  data  used  for  this  study  was  collected  by  two  magnetic 
spectrometers  aboard  Air  Force  Office  of  Aerospace  Research  satellite 
1969-25C  designated  OV 1-19  and  provided  by  the  instrument  builder  and 
prime  investigator,  A.  L.  Vampola,  Aerospace  Corp.,  El  Segundo,  Ca. 
The  purpose  of  this  study  was  to  find  the  frequency  of  occurrence 
of  REP  events  and  determine  the  suitability  of  this  type  of  data  for 
future  study. 

The  0V1-19  satellite  was  launched  on  18  March  1969  into  an 
elliptical  orbit  with  apogee/perigee  at  5790/467  km.  Inclination  was 
104.8°  and  the  period  154  min.  The  satellite  was  spin  stabilized  at 
8.4  rpm.  A  three-axis  magnetometer  tabulated  aspect  information. 

Four  hours  of  continuous  data  was  collected  and  transmitted  each 
tape  recorder  dump.  Over  the  life  of  the  satellite,  approximately 
25"  data  coverage  was  achieved. 

The  two  spectrometers  were  of  the  180°  magnetic  focusing  type 
similar  to  those  onboard  AF/OAR  satellite  0V3-3  and  schematically 
presented  in  Fig.  A-l .  (Vampola,  1969).  One  unit,  the  Low  Energy 
Magnetic  Spectrometer  (LEMS),  had  eight  electron  channels  from  53-444 
keV  with  geometric  factors  as  shown  in  Table  A-l  and  the  other,  the 
High  Energy  Magnetic  Spectrometer  (HEMS),  had  16  electron  detectors 
from  537  keV  to  5.1  MeV  as  shown  in  Table  A-2.  The  look  direction 
of  each  instrument  was  perpendicular  to  the  satellite  spin  axis  (for 
pitch  angle  scanning)  which  was  nominally  perpendicular  to  the  orbit 
plane.  Flux  counting  rates  were  sampled  once  per  second.  As  with 
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0V3-3  180°  Magnetic  Focussing  Spectrometer  (shown  without 
baffling)  from  Vampola,  1969. 
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TABLE  A- 


Low  Energy  Magnetic  Spectrometer  Channels,  Energy,  and 
Geometric  Factors. 
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TABLE  A- 


High  Energy  Magnetic  Spectrometer  Channels,  Energy, 
and  Geometric  Factors. 
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the  0V3-3  units,  the  flux  determinations  were  limited  at  low  counting 
rates  by  statistical  fluctuations  and  at  high  counting  rates  by 
telemetry  quantization  (Vampola ,  1969).  In  addition,  a  rectangular 
energy  versus  geometric  factor  response  is  assigned  for  each  actual 
channel  response  for  data  analysis  as  illustrated  in  Fig.  A-2. 
Instrument  response  to  an  input  spectrum  zero  everywhere  in  the  loss 
cone  and  isotropic  at  some  level  outside  the  loss  cone  is  shown  for 
Channel  4  of  the  low  energy  unit  in  Fig.  A-3.  An  order  of  magnitude 
drop  in  flux  response  roughly  occurs  7°  into  the  loss  cone. 

A  typical  high  to  low  altitude  pass  for  Channel  7  of  the  low 
energy  unit  is  shown  in  Fig.  A-4.  Inner  zone  (A),  slot  region  (B), 
outer  zone  (C),  and  polar  reqion  (D)  are  distinctive  features  of  the 
pass.  The  large  variation  in  count  rate  is  due  to  pitch  angle 
sampling. 

Fig.  A-5  for  an  0V3-3  sample  shows  the  outer  zone  at  low 
altitude  from  L=3-6  with  envelopes  drawn  for  the  top  and  bottom  of 
the  sampling  scatter.  The  top  envelope  is  representative  of  the 
locally  trapped  flux  and  the  bottom  is  the  flux  furthest  into  the 
loss  cone.  When  the  two  curves  touch,  as  in  Fig.  A-5  at  L=5.6,  this 
indicates  an  isotropic  pitch-angle  distribution,  as  predicted  for 
strong  pitch  angle  diffusion.  The  analysis  mode  is  then  to  view  all 
of  the  available  passes  for  this  envelope  closure.  Final  deter¬ 
mination  is  made  when  flux  vs  pitch  angle  plots  show  isotropic 
fluxes  filling  the  upward  looking  (downward  flowing)  loss  cone. 

In  order  to  accomplish  the  task  of  analyzing  the  data  set 
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FIG  A-2.  Geometric  factor-energy  response  function  assigned  each 
channel  for  data  analysis.  Total  area  outside  of  the 
rectangle  is  about  15%  but  both  functions  contain  equal 
areas.  A  slight  bias  towards  higher  energies  for  very 
steep  spectra  is  produced  (Vampola ,  1969). 
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FIG  A-3.  Instrument  response-pitch  angle  output. 
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FIG  A-4.  Analog  record  of  376+28  keV  electron  data  plotted  as 

counting  rate  versus  time  on  a  pseudo-logarithmic  scale. 
Universal  time,  magnetic  field  intensity,  invariant 
latitude,  L,  altitude,  and  MIT  are  listed  for  200-sec 
intervals.  Pitch-angle  distribution  sampling  produces 
the  large  variations  in  counting  rate  in  several  regions 
of  the  plot  (Vampola ,  1971). 
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FIG  A-5.  Plots  of  instantaneous  flux  measurements  versus  L  for  the 
712±137  keV  electron  channel  from  the  OV3-3  magnetic 
spectrometer.  The  scatter  in  the  data  points  is  due  to 
pitch-angle  sampling.  Envelope  curves  at  the  maximums 
and  minimums  are  shown  to  emphasize  the  data  that  are 
isotropic  in  pitch  angle  (Koons  et  al ,  1972). 
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as  objectively  as  possible,  a  set  of  criteria  was  developed  for  REP 
event  identification.  These  are  illustrated  in  flow  chart  form  in 
Fig.  A-6.  The  first  criteria  was  that  only  data  in  the  L  range  4-7 
would  be  studied.  A  total  of  320  passes,  clustered  as  high  as  4 
passes  on  a  given  day,  was  extracted  from  a  data  set  spanning  from 
day  77  of  1969  to  day  13  of  1970  (301  days).  Because  of  the  finite 
width  of  the  field-of-view  (-10°)  and  the  fact  that  the  satellite 
orbit  precluded  pitch  angle  sampling  to  0°  (and  180°)  everywhere  in 
the  orbit,  it  was  determined  that  data  below  1000  km  would  yield  loss 
cones  large  enough  for  a  high  frequency  of  complete  loss  cone  sampling. 
This  criteria  eliminated  88  cases.  The  third  criteria  demanded  that 
REP  fluxes  be  high  at  higher  energies  (1  MeV  ±  500  keV).  Seventy-two 
cases  failed  on  this  criteria.  Next  we  determined  that  if  we  had  data 
at  least  10°  into  the  loss  cone,  we  could  identify  real  precipitation 
based  on  the  channel  response  previously  described.  This  eliminated 
39  cases  of  which  17  were  less  than  10°  into  the  loss  cone  and  22  did 
not  sample  the  loss  cone  at  all.  In  some  cases  (19),  count  rate 
sampling  and  pitch  angle  modulation  were  phased  in  such  a  way  as  to 
make  event  evaluation  impossible.  These  cases  were  labeled 
"instrument  problems"  and  eliminated.  Finally,  103  cases  remained 
in  which  an  analysis  could  reasonably  be  made.  These  passes  are 
shown  as  passes  per  day  versus  day  of  the  year  in  Fiq.  A-7.  Of  these 
103  "clean  cases",  only  three  exhibited  Isotropic  fluxes  into  the 
loss  cone  at  relativistic  energies  (Identified  by  arrows  on  Fig.  A-7). 

The  REP  of  14  May  1969  (day  134)  Is  of  particular  Interest 
since  it  exhibited  the  "hardest"  spectrum  of  the  three  REPs  (spectra 
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FI6  A-6.  Flow  chart  Illustrating  process  of  elimination  through  the 
0V1-19  data  set. 
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FIG  A-7.  Outer  zone  passes  per  day  versus  day  of  the  year  for  the 
0V1-19  data  set.  Arrows  Indicate  the  three  event  days. 
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will  be  compared  later)  and  falls  in  a  period  where  Alaska  forward 
scatter  data  exists.  Figures  A-8,  a  and  b,  show  plots  of  differential 
electron  fluxes  versus  U.T.  (in  seconds)  for  the  various  energy 
channels,  LEMS  Bremsstrahlung  (background),  and  LEMS  protons  (>55  keV) 
for  day  134.  The  higher  energy  electron  channels  exhibited  low  fluxes 
and  are  not  shown.  The  outer  edge  of  the  outer  zones  at  local  morning 
(07:20  MLT:  79642-79659  U.T.)  and  local  midnight  (00:18  MLT ;  80184- 
80225  U.T.)  exhibit  envelope  closure.  The  local  morning  outer  zone 
exhibited  very  high  trapped  fluxes  (see  Fig.  A-9)  at  all  energies  but 
the  loss  cone  (-55°)  was  not  sampled. 

The  pitch  angle  distribution  for  the  local  midnight  (80196- 
80225  U.T.)  outer  zone  transit  are  shown  in  Fig.  A-10  for  energy 
channels  192,  537,  822  keV  .  The  nominal  loss  cone  is  ~66°  and  the 
symbols  "+/ou  represent  instrument  viewing  positions  up/down  the  local 
magnetic  field  line.  The  plots  show  high  fluxes  in  the  upward  viewing 
loss  cone  at  several  instances  in  all  energy  channels  while  the  down¬ 
ward  viewing  loss  cone  generally  show  normal  trail-off  to  lower  values. 
The  large  variation  in  precipitating  fluxes  is  due  either  to  the 
transition  of  the  satellite  through  regions  of  strong  and  not  so 
strong  precipitation  or  spatial  variations  in  trapped  flux  during  the 
30  seconds  duration  of  the  pitch  and  angle  plot.  The  plots  represent 
approximately  seven  complete  (0-180°)  pitch  angle  distributions  each 
and  537  keV  channel  has  each  data  sample  labeled  with  the  last  three 
digits  of  the  U.T.  in  seconds.  Deep  in  the  loss  cone  two  orders  of 
magnitude  variation  occur  between  210  and  218  whereas  a  half  order 
of  magnitude  between  205  and  220  or  201  and  197  occur  in  the  trapped 
distribution. 
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FIG  A-8a.  Differential  electron  fluxes  versus  UT  (sec),  B  (Gauss), 
L,  LAT  (deg),  ALT  (km),  for  LEMS  Proton,  LEMS  BREM,  822, 
537,  444,  376  keV  Electrons  for  DAY  134. 
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FIG  A-8B.  Differential  electron  fluxes  versus  UT  (sec),  B 
L,  LAT  (deg),  ALT  (km),  for  312,  250,  192,  139, 
53  keV  Electrons  for  DAY  134. 
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FIG  A-9.  Electron  differential  energy  fluxes  versus  UT  (hr),  LT  (hr), 
ALT  (km),  L,  for  192,  537,  822  keV  electrons  for  DAY  134 
(14  MAY  1969). 
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FIG  A-10.  Flux  versus  Pitch  Angle  for  electron  channels  192,  537, 
and  822  keV  for  DAY  134.  Flux  upper  limit  in  channel 
822  keV  is  2.696  x  101  since  the  high  count  rate  meter 
coax  cable  from  the  HEMS  to  the  satellite  data  system 
broke  during  launch. 
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A  time  sequence  in  pitch  angle  versus  differential  flux  for 
the  same  local  midnight  outer  zone  pass  (day  134)  for  channels  444, 
and  537  keV  are  shown  in  Figs.  A-ll  a  and  b  respectively.  Data 
samples  are  subjectively  connected  to  display  the  most  "normal" 
distribution.  This  is  difficult  because  of  the  wide  spread  in  pitch 
angles  (-50°)  between  samples  and  compounded  by  a  few  missing  samples. 
The  444  keV  channel  is  the  highest  energy  measured  by  the  LEMS  whereas 
the  537  keV  channel  is  the  lowest  for  the  HEMS.  Thus,  Fig.  A-ll  a  is 
165°  out  of  phase  with  Figs.  A-llb.  Each  sample  is  identified  by 
the  last  three  digits  of  the  U.T.  in  seconds.  "Normal"  upward  and 
downward  loss  cone  distributions  are  evident  at  all  energies  after 
sample  213.  The  upper  envelope  (dashed  line  -U)  connects  the  locally 
trapped  flux  (pitch  angles  -90°)  and  the  lower  (dashed  line-L)  connects 
the  lowest  flux  into  the  loss  cones  (pitch  angles  ~0°,  180°).  It  is 
evident  that  the  envelopes  nearly  close  in  at  least  a  few  places 
towards  the  outer  edge  of  the  outer  zone  though  the  instrument  is 
viewing  almost  2  1/2  orders  of  magnitude  into  the  loss  cone.  Most 
notable  is  the  isotropic  (upward  viewing  loss  cone)  flux  in  Fig.  A-llb 
at  sample  188  (80188  U.T.). 

The  magnetic  activity  for  this  event  as  well  as  Alaska  forward 
scatter  absorption  are  illustrated  in  Fig.  A-12.  The  Alaska  forward 
scatter  absorption  on  day  134  was  greater  than  10  db  from  19:40  to 
02:10  U.T.  (09:40-16:10  Alaskan  local  time)  with  a  peak  absorption 
of  19  db  at  20:12  U.T.  The  only  other  absorption  recorded 

was  early  on  days  132  and  137.  Dst  shows  disturbed  conditions 
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FIG  A-lIa.  Time  sequence  of  differential  electron  flux  versus  pitch 
angle  for  the  537  key  channel  on  DAY  134. 
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FIG  A-llb.  Time  sequence  of  differential  electron  flux  versus  pitch 
angle  for  the  444  keV  channel  on  DAY  134. 
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FIG  A-12.  Dst  for  days  132  through  137;  AE  Index  for  parts  of  the 
UT  day  for  days  134  and  135;  Alaskan  forward  scatter 
absorption  events  are  represented  by  the  hatched  boxes 
(Bailey,  Personal  Communication,  1973). 
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from  late  on  day  132  through  recovery  on  day  137.  A  positive  excursion 
is  particularly  noticeable  near  22  U.T.  on  day  134  near  the  time  of 
the  satellite  REP  observation.  This  is  followed  by  a  >100y  depression 
early  on  day  135.  AE  shows  a  1500y  increase  at  20  U.T.  on  day  134 
two  hours  before  the  satellite  observation.  It  is  interesting  to 
note  that  although  day  133  and  135  are  strongly  disturbed  no  forward 
scatter  absorption  occurred  above  Alaska  (Thorne  and  Larsen,  1976). 

A  local  early  morning  and  late  evening  pass  by  0V1-19  occurred  on 
day  135.  The  local  early  morning  pass  did  not  sample  the  loss  cone 
but  even  though  the  evening  pass  did,  no  REP  was  observed. 

All  three  REPs  observed  occurred  between  local  midnight  and 
local  dawn.  The  expected  frequency  of  occurrence  of  REPs  based  on 
this  study  is  1-13%  with  a  90%  confidence  interval.  The  50%  con¬ 
fidence  level  indicates  a  fequency  of  5%  (see  Fig.  A-13).  This  is 
consistent  with  5%  frequency  for  moderate  (>8  db  absorption)  REPs 
found  by  Bailey  (1968). 

The  energy  spectra  for  the  three  0V 1-19  REPs  are  compared 
with  the  0V3-3  REP  measured  by  Vampola  (1971)  in  Fig.  A-14.  It  is 
evident  that  they  are  comparable  in  hardness  (approach  to  isoenergetic) . 
So  it  is  not  surprising  that  the  REPs  will  have  a  similarly  predicted 
effect  on  the  deposition  of  energy  (represented  by  ion  production) 
in  the  atmosphere  as  shown  in  Fig.  A-15. 

The  conclusions  of  this  study  are  that  the  OVI-19  was  capable 
of  observing  strong  diffusion  events  (REPs)  from  a  few  hundred  keV 
to  1  keV  energies.  More  REPs  would  probably  have  been  observed  if 
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the  loss  cones  were  sampled  more  frequently.  A  similar  data  set  with 
better  statistics  would  yield  more  morphological  information  and 
possibly  a  clue  to  the  type  and  relative  importance  of  various 
mechanisms  for  causing  REPs. 
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FIG  A-13,  REP  frequency  of  occurrence  for  the  0V1-19  Data  Base. 
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FIG  A-14.  REP  energy  spectra  for  the  three  0V1-19  events  in  1969 
compared  with  the  Vampola 0V3-3  spectrum  of  1971. 
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FIG  A-15.  Ion  production  rate  versus  altitude  for  the  0V1-19  REP, 

solar  radiation,  galactic  cosmic  rays  and  the  solar  proton 
event  of  1972. 
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Appendix  B  -  S3-3  Data  Collection 

The  primary  data  set  used  in  this  dissertation  was  collected 
by  a  magnetic  spectrometer  aboard  Air  Force  Space  Test  Program 
satellite  STP  74-2  designated  S3-3  launched  in  the  suniner  of  1976. 

The  S3-3  orbit  is  highly  elliptical  (schematically  shown  in  Fig.  B-l) 
with  apogee/perigee  nominally  8000/250  km.  The  maximum  apogee  of 
8035  km  was  reached  on  3  August  1976  and  continues  to  decrease  at  a 
rate  of  0.35  km/day.  Orbit  inclination  is  97°  with  a  period  of 
180  minutes.  The  nodal  drift  rate  is  0.27°/day  westward,  i.e.,  the 
local  time  of  orbit  nodal  crossing  decreases  0.049  hr/day  from 
launch  values  of  10:30  (ascending  node)  and  22:30  (descending  node) 
as  shown  in  Fig.  B-2,  panel  a.  The  line  of  apsides  (major  axis  of 
orbit)  was  initially  placed  less  than  90°  from  equator  on  the  sun 
side  (Fig.  B-l)  and  slowly  rotates  toward  the  dark  side  at  the  rate 
of  0.94°/day  from  a  launch  value  of  63°  as  shown  in  Fig.  B-2,  panel  b. 
The  satellite  spin  rate  was  maintained  at  approximately  3  rpm  with 
the  spin  vector  perpendicular  to  the  orbit  plane.  In  order  to 
maintain  the  solar  cells  facing  toward  the  sun  (and  the  backside  away 
for  thermal  control),  the  satellite  was  yawed  180°  on  days  40-48 
(around  14  February)  1977. 

The  S3-3  scientific  payloads  are  listed  with  brief  descriptions 
in  Table  B-l.  Of  prime  interest  is  the  CRLS-217  which  will  be 
described  here  in  detail.  Data  from  0NR-104  and  CRLS-218  are 
referred  to  in  Chapter  2  and  only  brief  details  will  be  provided. 
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S3-3  Orbital  parameters. 
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FIG  B-2.  Day  of  the  year  versus  a)  sun  time  of  the  orbit  phase 
and  b)  latitude  of  apogee. 
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TABLE  B-1.  S3-3  Scientific  payloads  with  prime  investigators 
in  parenthesis. 


The  satellite  payload  configuration  is  shown  in  Fig.  B-3  with  the 
spin  vector  perpendicular  to  the  page.  TheCRLS-217  location  is 
such  that  penetration  thresholds  are  approximately  10  MeV  for 
electrons  and  100  MeV  for  protons. 

CRLS-217  is  a  twelve  differential  energy  channel  magnetic 
electron  spectrometer  with  an  eight  channel  integral  proton  and 
alpha  (1  channel)  telescope.  Permanent  magnets  (.53  and  1.64 
kiloGauss)  momentum  analyze  electrons  entering  the  instrument  aperture 
focusing  them  onto  12  silicon  semiconductor  detectors  at  the  primary 
focus  (180°).  The  instrument  configuration  is  similar  to  the  0V3-3 
and  OVl-19  instruments  described  in  Appendix  A  except  that  the  Low 
Energy  Magnetic  Spectrometer  (LEMS)  and  High  Energy  Magnetic 
Spectrometer  (HEMS)  are  built  into  one  unit  as  shown  in  Fig.  B-4. 
Pulseheight  analysis  of  the  energy  deposits  in  the  detectors  is  used 
to  discriminate  against  bremsstrahlung,  noise,  and  penetrating  ion 
background.  Also,  an  additional  shiedled  detector  is  used  to  monitor 
the  background.  The  electron  output  of  each  detector  channel 
corresponds  to  the  uni-directional  electron  intensity  of  a  given  energy 
bin  (typically  ±  15  keV  for  LEMS,  ±  100  keV  for  HEMS)  and  from  these 
outputs  an  energy  spectrum  can  be  obtained. 

The  data  format  is  presented  in  Table  B-2.  Included  are  the 
identification  code,  channel,  sample  rate,  bits/sample,  and  geometric 
energy  factor  for  28  CRLS-217  and  2  (21,  22)  CRLS-218  measurements. 
Measurements  1-5  are  the  5  LEMS  channels  while  6-12  are  the  HEMS, 
with  all  channel  accumulators  read  out  every  1/16  second.  Note  that 


FIG  B-3.  S3-3  Payload  Configuration. 
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FIG  B-4.  S3-3  HEMS/LEMS  Layout. 
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TABLE  B-2.  S3-3  HEMS/LEMS  Data  Format. 
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the  geometric  energy  factor  decreases  with  energy  for  the  LEMS 
and  HEMS,  i.e.,  sensitivity  decreases  with  increase  in  energy  (given 
flux  will  produce  less  counts/sec  for  a  reduced  geometric  energy 
factor).  In  addition,  S3-3  HEMS  channels  are  more  than  a  factor  of 
2.5  higher  in  geometric  energy  factor  than  0V1-19  LEMS  channels 
to  ~450  keV  (see  Table  A-l)  but  more  than  a  factor  of  5.8  lower  than 
OV 1-19  HEMS  channels  to  ~1.1  MEV  (see  Table  A-2).  Measurements  13-19 
and  20  show  specs  for  the  integral  proton  and  alpha  channels  which 
were  also  sampled  at  a  rate  of  1/sec. 

A  functional  block  diagram  of  the  instrument  is  shown  in 
Fig.  B-5a.  The  lower  level  discriminator  is  set  at  50%  of  the  minimum 
expected  energy  deposit  to  reject  bremsstrahlung  and  noise  and  the 
upper  level  discriminator  is  set  at  an  upper  threshold  of  110%  of 
the  maximum  expected  energy  to  reject  penetrating  particles  (cosmic 
rays  and  energetic  protons). 

Each  electron  energy  channel  is  calibrated  to  find  the  channel 
response  curve  from  which  the  absolute  peak  (EAPEAK),  minimum  (EAMIN), 
and  maximum  (EAMAX)  values  are  determined  (see  Fig.  B-5  b).  A 
triangle  which  best  matches  the  slopes  of  the  response  curve  then 
yields  its  peak  (E),  minimum  (ERMIN),  and  maximum  (ERMAX)  values. 
Center  values  of  the  two  triangle  slopes  determine  EMIN  and  EMAX. 
Vertical  lines  are  estended  through  EMIN  and  EMAX.  The  vertical 
extension  of  these  two  lines  continues  until  the  resulting  rectangle 
area  formed  with  line  segment  H  equals  the  area  under  the  actual 
response  curve.  This  representative  rectangle  is  then  used  as  the 
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FIG  B-5.  a)  Electron  Channel  component  schematic j  b)  general 
response  curve;  c)  channel  values  for  E°  with  100" 
efficiency. 
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channel  response  for  data  analysis.  (This  is  identical  to  the 
description  in  Appendix  A  and  depiction  in  Fig.  A-2  for  instruments 
0V3-3  or  OVi-19.)  Values  for  the  above  parameters  for  an  assumed 
spectrum  of  E°  with  1 00%  efficiency  are  listed  in  Fig.  B-5  c.  A 
typical  three  dimensional  response  is  shown  in  Fig.  B-6  for  the  435  KEV 
channel . 

The  instrument  loss  cone  response  for  655  keV  channel  is  shown 
in  Figs.  B-7  a-c  as  the  result  of  an  input  distribution  isotropic 
everywhere  outside  the  loss  cone  and  falling  off  about  an  order  of 
magnitude  per  2°  in  the  loss  cone  for  pitch  angles  of  8°,  18°,  28°, 
respectively.  Response  to  a  sharp  cutoff  can  be  demonstrated  by  a 
zero  input  distribution  at  the  loss  cone  but  is  not  represented  here. 
Except  for  pitch  angles  less  than  -10°,  an  order  of  magnitude  drop 
in  response  is  about  9°  into  the  loss  cone.  This  response  obviously 
affects  the  subjective  selection  of  isotropic  distributions  in  the 
loss  cone  and  is  further  discussed  in  Appendix  C. 

The  spin  rate  of  -3  RPM  allows  the  instrument  to  complete  a 
pitch  angle  plot  (0-180°)  in  approximately  10  seconds  (compared  to 
~4  seconds  for  0V 1-19).  At  16  samples/sec,  the  instrument  will  collect 
-150  samples  for  a  pitch  angle  plot,  nearly  a  factor  of  40  better 
than  the  0V 1-1 9 .  The  S3-3  advantage  is  clear  when  the  two  data  sets 
are  compared  for  analysis. 

The  proton/alpha  telescope  consists  of  two  silicon  detectors 
in  series  as  shown  in  Fig.  B-8.  Detector  D1  (100  pm  thick)  has  a 
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FIG  B-6.  Typical  three  dimensional  response  curve  for  HEMS  435  keV 
electron  channel. 
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FIG  B-7.  Typical  instrument  loss  cone  response  for  the  655  keV 
electron  channel  for  pitch  angles  8°,  18°,  and  28°. 
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lower  threshold  of  50  keV  and  detector  D2  (200  urn  thick)  has  a  lower 
threshold  of  100  keV  and  is  the  anti  to  reject  energetic  protons  and 
cosmic  rays  penetrating  through  01.  Particles  depositing  more  than 
about  4  MeV  are  assumed  to  be  alphas  since  p+  >3.2  MeV  would  pass 
through  D1  and  be  antied  out  by  02.  All  channels  are  sampled  at 
4  samples/sec  and  E  80  keV  and  E,>  4.0  MeV  are  also  sampled  at 
1/sec. 

The  S3-3  wave  experiment  by  the  Berkeley  group  (Mozer  et  al . , 
1977)  obtained  the  three  orthogonal  components  of  the  static  or 
fluctuating  electric  field  by  measuring  the  potential  differences 
between  three  orthogonal  pairs  of  separated  spheres.  Two  pairs  of 
spheres  were  separated  by  37  m  in  the  satellite  spin  plane  and  one 
pair  by  6  m  along  the  spin  axis.  One  sphere  on  the  long  boom 
could  be  operated  as  a  Langmuir  probe  to  make  relative  density 
fluctuations.  Broadband  signals  (.05  -  16  kHz)  were  transmitted  in 
real  time  to  ground  stations  for  15  min  periods  on  about  25%  of  the 
orbits  (Kintner  et  al . ,  1978).  The  percentage  of  concurrent  real  time 
wave  data  and  stored  plasma  data  (CRLS  217)  is  much  less.  Only 
filtered  wave  data  in  bands  from  .03-100  kHz  are  available  for  most 
of  the  stored  plasma  data  presented  here. 
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S3-3  Proton/alpha  telescope  data  schematic. 
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Appendix  C  -  REP  Event  Identification  and  Catalog  Use 

In  order  to  select  REP  events,  it  was  necessary  to  be  able  to 
identify  regions  of  isotropic  fluxes  on  microfiche  time  plots  of  the 
electron  and  proton  fluxes  indicative  of  strong  diffusion.  Fig.  C-l 
shows  such  a  plot  for  the  12,  235,  435,  655  and  850  keV  electron  and 
80  keV  proton  channels  for  rev  1661  on  a  logarithmic  scale;  UT  is 
measured  in  sec.  The  microfiche  displays  only  every  fourth  data  point. 
The  variation  in  fluxes  are  due  to  pitch  angle  modulation.  The  key 
feature  one  must  recognize  is  the  missing  upward  viewing  loss  cone 
indicative  of  a  full  or  isotropic  loss  cone.  This  feature  is  readily 
evident  throughout  the  12  keV  electron  channel  from  L  =  4.0  -  8.0  and 
the  80  keV  proton  channel  from  L  =  4.5  -  6.0.  The  235,  435,  655,  and 
850  keV  channels  exhibit  missing  loss  cones  at  L  ~  4.5  with  the  rest 
of  the  outer  zone  mostly  erroded.  Effects  such  as  data  drops  and  low 
altitude  diminished  flux  values  similarly  appear  as  missing  loss  cones 
so  it  is  extremely  important  to  follow  up  with  examination  of  the 
actual  pitch-angle  plot.  Fig.  C-2  shows  pitch-angle  plots  for  the 
235,  435,  655,  and  850  keV  channels  of  the  example  in  Fig.  C-l  for  UT 
61092  and  61109  secs.  The  magnetic  field  direction  is  zero  in  the 
southward  direction  at  the  equator  so  for  this  Southern  Hemisphere 
example,  180°  is  the  up-the-field-line  direction.  The  pitch  angle  is 
approximately  20°.  Also,  every  data  point  is  displayed  and  differ¬ 
ential  flux  values  can  be  obtained  by  dividing  by  the  geometric  value 
per  channel  (Table  B-2).  The  top  panel  shows  normal  loss  cone  dis¬ 
tributions  both  up  and  down  the  field  line  indicative  of  expected 
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radiation  belt  distributions.  (However,  the  belts  are  disturbed 
since  normally  quiet  times  produce  loss  cones  3  1/2  orders  of 
magnitude  down  from  trapped  flux).  The  bottom  panel  (data  toward  the 
outer  edge  of  the  outer  zone)  shows  the  pitch  angle  17  secs  later. 
Although  normal  distributions  are  evident  down  the  field  line, 
viewing  up  the  field  line  shows  isotropic  distributions  in  all  channels. 
Similar  plots  can  be  obtained  for  the  proton  channels. 

Since  the  satellite  spin  axis  is  not  always  directly  per¬ 
pendicular  to  the  local  magnetic  field  line  and  the  instrument  has  a 
field  of  view  of  approximately  5°  (half  angle),  we  chose  only  those 
cases  with  data  at  least  10°  into  the  local  loss  cone.  A  few  exceptions 
were  made  when  the  excursion  into  the  loss  cone  was  greater  than  6° 
and  isotropic  fluxes  were  me-sured  on  previous  and  succeeding  revs. 

Less  than  3%  of  the  325  events  selected  for  further  pitch-angle 
analysis  were  rejected  for  the  'less  than  10°  in  the  loss  cone' 
criteria. 

Apogee/perigee  bias  in  the  flux  data  was  investigated  and 
found  to  be  not  important  except  in  the  actual  appearance  of  the 
flux  versus  time  distributions.  At  apogee  the  loss  cones  are  small 
and  the  counting  rates  are  high.  Here  we  merely  insisted  on  the 
'10°  into  the  loss  cone'  criteria  for  validation.  At  perigee  the 
loss  cones  are  large,  the  extent  (in  time)  of  the  entire  outer  zone 
is  narrow,  and  the  flux  counting  rates  of  both  trapped  and  precipitated 
particles  is  small.  The  learning  process  to  identify  isotropic  fluxes 
on  the  microfiche  plots  was  accomplished  by  a  series  of  trial  and  error 
analyses  taking  about  3  months  of  time  and  150  samples.  A  high 
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Fig.  C-l.  Logarithmic  flux  (counts/sec)  versus  UT  (sec)  for 
the  12,  235,  435,  655,  850  keV  electron  and  80  keV 
proton  channel  for  rev  1661. 
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Fig.  C-2.  Flux  (counts/sec)  versus  pitch  angle  (degrees)  for 
rev  1661  at  61092  UT  (top)  and  61109  UT  (bottom) 
for  the  235,  435,  655,  880  keV  channels.  Loss  cone 
is  ~20°  with  the  upward  looking  loss  cone  at  160° 
indicated  by  a  vertical  line. 
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Fig.  C-3.  Microfiche  layout  with  each  logarithm  of  flux 
(counts/sec)  as  ordinate  for  each  channel  and 
abscissa  the  UT  (sec)  every  180  sec  for  900  sec 
(16  samples/sec),  ALT  [km],  LONG  [deg-E], 

LAT  [deg],  B  [Gauss],  L  [  ],  ILAT  [deg],  MLT  [hr]. 
E.G.  Eight  panels  (15  min  each)  of  12,  33,  70 
electrons,  8  panels  110,  160,  235  keV  electrons, 
etc.  P  =  Protons,  A  =  Alphas,  PD  =  Proton  Detector, 
BG  =  Background. 


232 


233 


Table  C-l.  Listing  of  individual  microfiche  by  Readout  REV  NO 
(Station  and  Number),  DATE  of  REV  (Last  Digit  of  Yr, 
Month,  Day),  Beginning  and  Ending  TIME  OF  REV  IN 
HRS  (Hours,  Minutes,  Seconds),  Beginning  and  Ending 
TIME  OF  REV  IN  SECS  (Seconds),  Precip  Category 
(1  <  235  keV,  2  <  435  keV ,  3  <  655  keV ,  4  <  850  keV ) , 
COMMENTS  (12  =  12  keV  electron  channel,  etc.;  80, 

150  =  80,  150  keV  Proton  Channel,  etc;  L  =  L  value). 
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DATE  TIME  Of  REV  TIME  Of  REV  I 

Of  REV  IN  HRS  IH  SECS  | _ _  _ _ 

.  *091*  9*02*  11*011  3482*  *0011  4  nt  Jf< '  C  *r-<*  *■  '  •"  A+4 

UW<  - HU  OW1 - 00*1  9  l 

i  Mil*  2 5151*  552*  05750  000*  \$3,liOfe  trfO 

I  0—00  UU — 01 100  0000  MU  . 

1*0920  554*  *1359  *0*4  0059  t  #*.,»*«- A  >k~  y  3  -  f  L  *  7  , 

!  *0920  75542  01151  *75**  *0**1  1  ll  JJC  £<*-*  *  »  ,  f 

.  *9920  1*09*0  »»»  000*9  —00 

1*0920  1*0921  15205*  *09*1  5*25*  £  /f.  9  »1  *  >*>’.»•  *  7*,  W,  /pR'  *  * 

i  *<Kt0  152110  >9074*  omo-*m*  „  _ 

.*0920  152115  170745  *5275  *1*45  ,V<’fc-'0  iO  -1 

»  40000  21102*  000000 — 51700  OHIO — 

.  $422p  211053  *50520  7C713  051*0  /2,  jjc  A  O- »»  Jlrf  Fj 


*£320  211053  *50520  7^715  051*0  Jl,  JJf  *  ■*  <r»®  i»f  (*<P  i*’!  ^  *  *V 

11017  7**0  HQ99 

60921  2042*  51010  74**  11090  J  £  »**-/•  O'  ff-.  >3?* 

♦0901  514»5 — 5004* — 1*915-  -tOlH — 1  7  ' 

*0921  5105*  5024*  11914  161*6  f  '7,33.*  i r  >*-  2 3  -  O-J**'/ ‘  *  *0  P  '  C* 


BOS  0602  *0921  110042  121509  2 
000-0*4*  *0971  12*00*  MWM  < 
BOS  0601  60921  121530  135902  * 
■OOP  0*00  *0051  *2**10  >010  I 

BOS  0609  $0921  225539  1552  < 

WO  4*90  *000*  *555 — *54*4 — 

XX>  0600  *0922  1555  15410 

COO  0*12  *092*  164611  175944  < 

BOS  0*11  *0922  100025  194*5*  < 

BOS  0*15  *0923  14426  50*16 

6UA  0614  *0925  50*5*  44252  1 

COO  0*19  *0925  104041  115220  2 

BOS  0617  *0925  115257  155907  * 

BOS  0*25  *0923  225545  254750  I 

CUA  0*22  *0925  25474*  15409  ( 

HUl  0*29  *0924  192700  203050  2 

COO  0620  *0924  203919  222 527  1 

COO  0*4 j  7792*  11454  22*07 

BOS  0639  *092*  22624  41259 

BOG  0*46  *0?27  10424  24921 

BOG  064*  *0927  2493*  40247  ] 

R06  0651  *0927  115404  134120  < 

HUl  0*49  60927  154137  153240  * 

CUA  0654  *0927  215515  251321  2 

COO  0653  *0927  231350  5350  « 

CUA  0*59  *79?0  135401  15374*  1 

HUl  0650  *0920  155001  173254  ! 

COO  0*65  60929  *410*  7515*  1 

COO  06*4  60929  75215  95915  \ 

CUA  0*71  *0929  170951  105333  t 

COO  0**0  *0929  105350  204347  1 
COO  0*7*  *0930  152*43  164442  | 
OOS  0*7$  *0930  1*4511  102452  I 
COO  0*04  *1001  114540  132043  < 
BOS  0*03  *1001  132059  151544  - 
BOS  0*07  *1002  1207  13*5* 

HUl  0*0*  *1002  13720  51029 

BOS  0690  *1002  90015  102557 

XtO  0*09  *1002  102*10  120*35 
HUl  0*93  *1002  16041*  191539  < 
COO  0*92  *1AJ>2  19155*  210239  . 
XND  0*97  *1003  50025  42232 

ROC  0*95  *1005  42249  *5043 
COO  0700  *1005  145500  1*065* 
BOS  0*99  *1005  1*0*55  17S329 
BOS  0700  *TW*  114550  125710 
HUl  070*  UAA*  125735  144421 
HUl  0717  *1005  03*50  9550* 

If®  0713  MJU15  95522  115512 
COO  0721  fT00*  52742  *4544 

•OS  0720  4100*  *4*03  02*05 
HUl  0730  *1007  21029  35515 
BOS  0729  *1007  33551  01*35 

CUA  0735  *1007  21Q042  200* 

HUl  0  754  OTTO  2023  20*59 
PC  074$  *1000  195903  211027 
HUl  0741  *1000  211042  223721 
HUl  0754  Cm*  74239  05323 
BOS  0753  41A10  05339  104057 

COO  07*4  *1011  132050  144*5* 
XfO  0763  *1011  144713  1*2713 
COO  0772  *1012  101921  11372* 
HUl  0770  *ifiif  113740  131740 
BOS  0777  *1013  11232  25035 
•UA  0775  OjffU  23051  41040 

OUA  0704  *1014  10135  22*17 
ROC  0703  *1014  22*35  35951 

COO  0793  «lH5  0025  2002* 
COO  0791  *1015  20042  >4037 

BOS  079*  *1015  144343  1*1302 
BOS  0795  *1015  1*1119  173952 
CUA  0799  *1015  214017  223313 
COO  0790  *1013  225629  5031 

XK>  0003  *nn*  93429  111915 
HUl  0002  *101*  111932  123245 
COO  0005  *101*  171159  104150 
BOS  0009  *101*  10421*  201541 
•UA  0006  *TB1 7  2724  13051 

HUl  0007  41017  150*2  21*52 

•UA  0012  *1017  122129  111502 
INC  0011  *1017  111532  1*1943 
•OS  001*  *1017  211710  222024 
•UA  MIS  *1017  222041  1023 
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MV  NO. 

0AT2 

TIME  I 

or  rev 

Tim  or  ttv 

07  02V  IN 

MRS 

IN  52C3 

•OS 

0819 

*1016 

65513 

102241 

32113 

37341 

-A  f  V  .© 

*0  L  ->  p 

DC 

0810 

UAU  102256 

120936 

37378 

43776 

t.  iJ  «  ’  < 

,  TOjSs  1  '  /*P 

503 

882* 

61019 

406 

12204 

246 

4926 

■UA 

0623 

41019 

12222 

30224 

4942 

10944 

505 

0626 

61019 

64615 

101523 

31575 

34923 

12  A-  ■**»  "  J 

DC 

•627 

41019 

101123 

115603 

34663 

43063 

503 

0532 

61019  235236 

10213 

65954 

3733 

.  *  >/  . 

/©- J  *  7 

COO 

0631 

61620 

10230 

25049 

3750 

10249 

lit) 

4  t 

DC 

0635 

61020 

63526 

100620 

30928 

14360 

12  t*  <*'f 

p 

COO 

0634 

61020 

100637 

114630 

34  397 

42390 

•03 

0049 

61020 

234102 

5220 

65262 

3146 

/2  **7*/ 

*o 

MUL 

0634 

61021 

5236 

23915 

3156 

Cll 

il, k*f-n 

503 

0646 

61021 

62407 

94602 

30247 

15262 

n  i-f’ 

COO 

0645 

94616 

113455 

35296 

41695  4 

?}/  f  t-s  7  >- 

•03 

8651 

61022 

60456 

95417 

29044 

54577 

_J»‘  < 

. 

COO 

0650 

61022 

936  32 

112310 

54592 

46990 

ft  L  -  t  J3 

coo 

0661 

*1523 

75437 

•2405 

25477 

33445 

t  c-  ff/.f 

•03 

0660 

61023 

92422 

111161 

33642 

46261 

i  L-t 

BOS 

0669 

61024 

104726 

120704 

36646 

43426 

DC 

5667 

6  m_* 

120722 

135719 

43442 

•0239  , 

•OS 

0677 

6TTTT5 

100039 

111709 

36039 

40429 

503 

0875 

61025 

111726 

125039 

*0646 

44239 

503 

0664 

*1026 

101131 

115246 

34691 

42766 

503 

0663 

61026 

115305 

133300 

42765 

46780  J 

COO 

0693 

61*17 

41624 

64213 

25104 

31333 

BOS 

0691 

HO  2  7 

64230 

102222 

31150 

J*s« 

*  L-  k 

•OS 

0697 

6W75 

11504 

22610 

4504 

•770  f  go  f  *  L  x  *  *  > 

MUL 

0696 

41026 

22627 

41314 

6767 

15194  I 

COO 

0902 

61026 

100232 

112631 

34152 

41311 

COO 

0901 

61026  112646 

130641 

41326 

47321  | 

!  J3  L't 

?v 

COO 

6911 

41029 

92052 

110255 

33452 

14775 

coo 

0910 

61029 

110312 

125631 

39792 

44591 

coo 

0916 

61030 

91313 

110029 

33193 

39429 

/?  3>c* 

coo 

0917 

41030 

110133 

124421 

39493 

45641 

'_/<?  7** 

BOS 

0924 

61531 

91443 

105639 

33263 

39519  ^ 

' i. »f t-  ..  * . 

i.  -  7  j 

70S 

0923 

61031 

103654 

123157 

39534 

4*117  ; 

j 

•J  r*a  r*  c 

QUA 

0926 

61031 

175354 

192731 

44434 

70051  j 

* j  a  --  L  r . 

'  ro  rl  t  »  •  3 

MUL 

0927 

61031 

192747 

212710 

70067 

77230  | 

fo  ;j  l-r--5 

MS- 

MU- 

-4&C1 

9MM- 

it»ii  j 

. 


If®  0951  01101  104909  121912  30949  44SS2 
MUL  0930  01101  101022  193233  05702  70353 
MUL  0934  01J01  193249  211424  70369  70404 
•OS  0940  "5X102  90020  102013  32100  37173 
POC  0939  01102  102030  120025  37590  43505 
M>  0940  01102  1  00335  1  92120  05015  09000 
MOL  0943  ftXlP2  192144  210130  09704  75090 
•03  0949  01103  115359  131130  42039  07490 
•03  0940  01103  131153  145205  07513  >3525 
BOS  0957  01103  234730  10525  *5050  5925 
DC  0955  01104  11542  24542  3442  5442 
DC  0902  01104  233451  4115  54091  2495 
MUL  0900  *TT05  0152  23253  2512  5175 
MUL  0909  fllOS  232211  3310  54131  1595 
COO  0900  01106  3332  22011  2012  5411 
COO  5974  01106  111541  122651  40541  44011 
COO  0973  01106  122707  101347  44027  51227 
2ND  0979  5TTTJ7  20745  32535  7*05  1  2335 
QUA  0970  01107  32551  50540  12351  10340 
COO  0904  01107  140110  151235  50470  54755 
COO  0902  0*157  151251  105917  54771  01157 
DO  0900  01100  75114  90907  20274  12947 
•05  0907  01100  90923  104914  32903  30934 
IND  0990  mi9  43932  55720  10772  21440 
XMD  5994  01109  55930  73733  21570  27453 
•Of  1005  «H0  102311  114054  37391  42054 
DC  1004  OHIO  114110  132113  42070  40073 
DC  1011  SITU  11452  *3904  4*12  0544 
MUL  1009  OUU  2*421  41251  9541  15171 
QUA  1010  01111  220309  231417  79309  03057 
6UA  1017  01111  231432  10109  03072  3669 
•03  1021  HTl2  91000  103503  33000  30105 
DC  1020  01112  103522  12*0*7  50122  40127 
MUL  1025  01112  214940  230052  75500  52552 
•UA  1024  01112  230109  4740  02009  2000 
COO  1030  5TT13  120127  140052  45007  50012 
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2067 

6336  1 

,  .  1 

•UA 
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14553 

33246 

4353 

12766  ^ 

;  13 1?+  2 

COO 

1036  41114 
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44665 
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|  <?,»<•  3 

BOS 

1637  61114 
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49639 

55325 

fc.-C 

•09 
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91421 

162729 

13361 
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)1.  S*i  ;  •  •, 

iro 

1044  41115 

102751 

121416 

37471 

44056 

izin.’t  . 

•UA 

1050  41115 

216924 

222787 

74146 

•6627 

n  .*  v|_ 

•UA 

16*9  41115 

222724 

722 

66644 

442  , 

1 1  ?  1  -  *  •  >  . 

60S 

1059  4lTl6 

235 SSO 

11133 

•4030 

4293 

| 

•UA 

1057  nri7 

11149 

25144 

4309 

16364 

l  'ti.  >4. 

COO 

1063  41117 

114454 

131046 

42414 

47446 

|  ,7 

*05 

1041  41117 

131165 

144452 

47445 

51012  1 

1 

COO 

1672  rmo 

5*429 

70051 

20169 

25251  ItL-h 

MUL 

1671  41116 

75167 

•3424 

25267 

30644 

•08 

io7#  rrm 

111912 

125020 

45752 

44220 

DC 

1076  *1119 

1258*5 

141707 

44235 

51427 

•U* 

1690  6TT20 

175158 

151247 

41310 

45547 

n  >  ^ 

MUL 

1666  41120 

151*64 

199947 

45564 

71957 

/jr.  >yt  -  6y 

COO 

1695  T1T21  165127 

120994 

19057 

45744 

COO 

169*  41121 

126922 

1*4917 

41762 

•9757 

COO 

1153  Writ 

14221 

*6441 

4141 

1**41 

ill  L  '  U  f 

•OS 

1101  41122 

*64*7 

44614 

11217 

16514 

U3  L  * 

•09 

1111  4II21 

1622*6 

114444 

*7*50 

4*464 
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ttv  NO.  DAT*  Tim  or  »cv  tp*  or  »*v 

OF  RtV  IN  WS  W  MCS 


•os  ms  tii?i  wm 

BOS  11U  *1124  769SZ 

•os  1117  bu*4  6*eo4 
•os  lira  ants  mtm 

FOG  1125  411*5  52654 

pc  ini  aim  mm 

•UA  1150  aim  15935 

pc  nsa  aura  laaaia 
coo  nsa  aim  uiiia 

■OS  1144  a 112 7  122551 
•OS  1145  41127  135451 
BOS  1151  61124  41101 
BOS  1150  aim  102210 
•UA  115a  a1l29  204 

KUL  1154  41129  10454 

FOG  1142  41124  204411 
PC  1140  41124  220025 
PC  1172  fTTSO  1754  04 
•OS  1140  41150  164451 
COO  1177  4*701  142504 
BOS  1174  41201  154717 
FOG  1165  C1702  111009 
•UA  1162  41202  125417 
I  HD  1166  m02  250244 
MU l  1164  41205  1404 

BOS  1191  41205  125556 
BOS  1140  41205  145521 
PC  1144  Wife  14411 
•UA  1145  41204  50400 
ZHD  1205  41204  225512 
«UA  1202  41204  254425 
PC  1212  41175  192020 
COO  lt04  fciiiS  t 05755 
FOG  1219  41204  140722 
CUA  1216  172515 

Z»C  1225  41207  175504 
BOS  1224  4I2Q7  165627 
BOS  1251  41206  64520 
PC  1229  41206  75420 

HUL  1234  41206  163545 
I*C  1233  41206  144437 
BOS  1259  6T7&4  62609 
PC  1237  41209  75654 
COO  1246  mic  11453 
ZHD  1244  41210  45247 
PC  1252  {ITU  140 
HUL  1251  61211  12403 
GUA  1259  61211  204901 
PC  1257  4ini  Z 15929 
•UA  1247  41712  173504 
MUL  1244  4121*  144415 
MUL  1275  41213  171953 
•OS  1273  41213  1BS054 
MUL  1263  41*14  146414 
COO  1260  61214  151713 
COO  1269 '{Tils  105234 
•OS  1247  61215  120340 
COO  1294  TIT 16  73900 

FOB  1295  41214  90324 
MUL  1306  mi  7  93457 
FOG  1343  61*17  112231 
BOS  1311  nnt  65444 
MUL  1310  61214  M422 
HUL  1314  41214  154425 
COO  1313  41214  171306 
FOG  1317  11W  5412 

HUL  1315  41219  13437 
MUL  1322  61*19  94431 
COO  13*1  41*19  114134 
MUL  1531  m*0  923*4 
PC  1330  41220  105139 
•OS  1334  (Till  •2606 
FOG  1135  62*21  103241 
MUL  1336*61*21  163727 
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•OS  2913  1A2A3  192124  212544  70104  771*6 

COO  2919  70704  91224  110443  33144  40103 
COO  2425  7170*  44701  14432  24428  '31472 
COO  2428  71705  1437*4  163*23  *2476  *9723 
COO  2934  70706  11724  11144*  29*44  34885 

•UA  2442  78707  4*31*  6*042  17*99  24442 

COO  2044  JB207  1441*2  16341*  *2912  *44*5 
KM  2049  70708  44121  62304  16081  22486 
•OS  24*1  71788  1130*7  122820  378*7  44400 
IX)  24*1  70708  132*45  150425  4834*  S426* 
•UA  24*4  71708  144*31  *11757  71011  74677 
806  24*4  71794  113*1  31050  4411  114*8 

COO  24*7  71719  300*0  61711  110*0  22051 

m  24*1  71714  7*446  01404  *6416  3*444 

COO  (4*9  7*719  110127  120033  36507  632*3 
HUL  2461  70709  19*610  17**96  *7*79  64416 
806  2462  -70781  11*241  21*007  6746*  7S007 
HUL  2463  71711  4*41  22107  2741  0447 

HUL  2444  70710  34220  *171*  13*40  140*3 

HOG  2444  70710  438*2  19411  299*2  30975 

805  2447  70710  4*044  11*2*0  3*444  41570 

COO  2461  70710  12392*  142923  43545  1216* 
NUL  *444  70710  1*300*  172*37  3*80*  427*7 
KM  1970  70710  115029  202230  44429  733*0 

•OS  2972  70711  MU  21538  1041  §131 

•OS  2973  71711  41142  S1214  14502  117*4 

•OS  2474  70711  419S0  S0844  224*1  19)26 

806  247*  7S711  01S2*  111522  *1*26  14422 

COO  2474  71711  121*18  1411*3  44U*  *1*33 
•OS  2977  71711  1*1422  14*821  14242  41101 
•UA  2971  71711  114S10  14SS03  47*14  71713 
COO  2474  71711  203415  222451  74171  M44S 

HUL  2411  TSni  *4719  44114  11034  16876 
HUL  til*  71712  14*148  141804  *1441  tt6#9 
80S  *441  71711  1*89*  111141  MSI*  16641 
•OS  *44*  71711  141626  141*43  *0784  *7*21 

•OS  *49*  70711  19*422  21*64*  714*2  79*8* 

M  *996  7*713  *13944  *319  Ml 6*  1144 

•OS  2440  9V7T4  44907  646*0  17347  MM* 

80*  *999  79714  •*§**  942*4  304SI  *4974 
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«v  m.  date  tine  or  nrv  time  or  mv 
or  BCV  IN  KBS  IN  IEM 

COO  3000  70714  11105*  121*25  4025*  45545 

HUl  3002  70714  1645*0  103233  403*0  447*5 

•UA  3003  70714  1*3140  212*04  70300  77344 

IK)  3004  7071.4  222012  2530  000*2  1530 

406  3005  70715  12445  52212  5005  1211* 

*05  3007  70715  71757  *1521  24277  »S21 

3000  70715  101445  121154  17005  41*14 
HU l  300*  70715  133045  150030  40445  54510 

HUl  1010  70715  140741  100505  50041  45105 

INC  3012  -7*735  220054  2  35010  7*254  042*4 

005  3014  70714  35350  55117  14030  21077 

BOS  3015  70714  45110  04750  24470  11470 

BOS  3014  70714  102103  114427  17243  42247 

NUl  3017  70714  124335  14405*  45015  5205* 

•UA  3010  70714  154407  173732  54747  43452 

IK)  3020  70714  21131*  233040  77***  #4440 

006  3021  70717  2*52  22714  17*2  *034 

BOS  3023  70717  42254  02021  22974  10021 

BOS  3024  70717  *1504  111450  14504  40410 

COO  3025  70717  124*20  141332  44140  51212 

NUL  3024  70717  155313  171005  571*3  41005 

BOS  3027  70717  100*13  200437  45353  723*7 

BOS  3020  7Q 71 7  21552*  230313  70*2*  *2**3 

BOS  3030  70710  25057  4541*  10737  1777* 

COO  3031  70710  40530  75255  21*10  20375 

COO  3032  70710  *1410  104*20  14450  1**40 

COO  3033  70710  114042  114400  42122  4*540 

BUA  3034  70710  144514  144214  13114  40154 

POO  3035  70710  174144  1*3*0*  41704  7074* 


BOS  3036  70710  211443  223544  74403  01144 
BOS  3030  70719  23123  4204*  *003  1412* 

BOS  303*  7071*  14715  72525  20015  24725 
*06  3040  7071*  02434  102157  10274  17117 

COO  3041  7071*  112556  131031  41156  47*11 

006  3042  7071*  143540  161506  52540  5*506 
406  3044  7571*  201052  220013  72652  7*4*3 

406  3046  79770  20405  40120  7445  14400 

•OS  3047  70720  51621  45753  10*01  25071 

HUL  3040  70720  03053  *5422  J11SJ  15442 

BOS  304*  70720  110405  125052  19*45  44252 

BUA  3050  70720  143*10  154722  52750  54042 
•OS  3052  70720  1*4254  214010  70*74  70010 

HUL  3053  70720  234030  11540  *5710  1740 

BOS  3055  f0721  43231  62*4*  14151  2110* 

IK)  3056  70721  71111  *2620  27071  11*60 

HUL  305*  £0721  13303*  151*21  4B63*  55163 
NUl  3065  70722  *424*  114012  14*4*  40012 

BOS  3071  79725  10410  92120  5010  12000 

COO  3074  70723  121X22  140036  41602  56*14 
COO  3002  7B2F*  114100  134032  42106  4*fl2 
NUl  3007  70725  22533  42257  6713  15777 

COO  30*3  Z£Z£5  151022  170744  54422  41644 
•OS  30*4  7Q72fe  4535*  05121  1743*  24601 

COO  3106  70727  72225  *1*45  24545  13505 
BOS  110*  70727  161151  100*12  56311  45152 
•UA  1117  25226  154043  174107  54443  63667 
HUl.  311*  7072*  11652  50521  11012  16121 

COO  3121  7072*  52034  72557  1*716  24757 

NUL  3122  70729  *3020  112115  14220  40071 

NUL  3121  7072*  122204  1*173*  44524  51459 
OUA  5124  7672*  151644  171400  55004  42046 
•UA  3125  7072*  161317  20103*  655*7  7243* 
XMD  3124  7072*  211141  230110  74301  62074 

BOS  3120  797To  302*5  45*5*  10*45  17*94 

605  912*  70730  41124  75424  22204  26564 

COO  3130  70730  *0*26  105253  12*44  1*171 

COO  3111  70710  111156  11**1*  42716  4*714 
•UA  1111  76710  175644  1*4211  64724  76*11 
DC  1114  70730  20411*  223636  7447*  61116 
NUl  9115  76730  231751  11505  65071  5705 

BOS  3137  UTsi  1103*  72750  1*01*  24078 
•OS  3110  70711  04400  102420  11540  17440 

NUL  111*  70711  115092  112054  41112  46054 
HUL  3140  70731  1454*0  1*1723  53606  56443 
HUL  3144  TJSSl  15213  14*51  4753  137*3 

HUl  3155  79*02  71707  *142*  24227  13264 

606  3140  70602  215*24  2S5422  7*144  64102 

BOS  3144  Wl  12342*  1*365*  4536*  >2734 
606  314*  7WW  32350  52110  12230  1*270 
HUL  1171  70004  100407  200325  45147  72205 
606  317*  MOOS  64624  102621  11704  17561 

HUL  317*  70065  162640  121621  17406  41621 
•OS  3161  70601  111200  165221  55*20  40741 

606  3163  70005  204721  224516  74041  61*14 

COO  3164  T9S9*  232*  1*142  140*  6102 

606  316*  76004  51713  73436  26213  27274 
BOS  3107  70006  05050  103105  11050  STMS 

COO  3200  70004  111010  13272*  41410  4044* 
BOS  110*  70004  14155*  1*2150  *2554  1*010 

COO  31*1  76064  261*36  221*10  71170  Mill 
BOS  1161  79537  21235  4**44  7*15  14*64 

•OS  31*4  700*7  52035  7M19  1*235  *»TJ 

HUl  31*5  76007  M523  1**23*  2*123  Mil* 
COO  31*4  70007  111052  125*67  46*52  44747 
KM  91*7  76667  142100  15353*  11440  17554 
•Of  31*6  7Q607  16541*  16515*  6667*  67*1* 
HUL  1243  TVffB*  72205  *1*2*  265*5  SM44 


NUL  1200  76000  2*04*2  M  7*4*2 
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»ev  mo.  dati  rim  or  *tv  Tim  or  *cv 

O*  *ev  IN  M*S  IN  «C1  L _ _ _ _ _ 

COO  UM  70811  44050  11050*  39330  144«S  ]  fO,  » L  '  w  *T 

HUl  3229  70011  121*10  1*0130  **170  5049#  jij.lsi*/ 

HUl  3230  70011  1S003S  1*9733  9*033  *1073  4  tJ  t  .lJ*  r*  . 

6LU  32  3*  70112  19473*  17314*  9**94  *439* 

COO  3244  70013  93434  71420  20044  2*34* 

003  3244  7P013  104407  202411  *7447  737*3 

HUl  3234  71fl4  101231  121444  17131  441*4  2  9  T  w  *  (,  JO, 

•UA  1233  70014  170600  102014  *13*0  **444  4  IT  l  -  /  po,  ,0  i  to  ? 

OOS  1234  Tftftl*  1 44  744  212443  71244  77*03  i  J|f(.  )  .  .  .  t  *  *  *  ’ 

COO  325  7  7*013  12014  21020  4*1*  4*20  ,>n:|»  f  O  L  '  + 

COO  3240  70015  45745  03414  25043  12034  -  *  **  £•■<*>-*  fO.  '  »*,* 

COO  3263  25S15  1*42  34  203443  *  7334  74303  ,  }1  Po  L--  + 

COO  3264  70016  62003  02313  23203  1*115 

HUl  3274  75017  124133  140522  43443  90722  90  -Vi  > * 

I*C  32*0  70017  100033  144497  *4091  7121?  i3-'"  L s  r»j  J-*  - 

COO  1201  70£J 7  203704  **5427  7*424  *24*7  JIT  7-  -  ' 

006  3203  7001*  2449*  44711  10140  17233  •  pe  ~  >  i' 

COO  3205  70010  74330  4411#  27*30  14070  ♦.JJ.t  •!*•*,  23.'  i.  -  •  /  frV  t  "  - 

COO  32*4  7001*  113402  113423  41442  4*403  £  *J  .  ?X  2* rvi/^'fc  7." 

90S  32*7  700*0  143527  1*1241  02327  045*3  t  2ir  r3  l  >  v 

HUL  22*4  7001*  202023  222320  71703  *0720  ^  l»r  /  =<*  *0  »  10 

*03  3241  7®1T4  22054  41014  *434  19444  J3i7<  ^ird»6Ut*«  '  TO  i 

COO  1244  70014  7145*  41211  2*04*  11113  2  .  S  i  t  -  7  6  7-2  L  «0  .  k**  <-**  - . 

COO  3247  144131  213410  70413  77450  ^  |J(?r 

HUL  3303  70020  10233*  122091  17*1*  44451  7I  >3*-"?  #0  ,r£  tv? 

COO  3310  f5S?l  69755  *3311  25075  SIU1  *  J .  3  3  l  >  Tt-i*/ ‘  O  Ifl  if? 

COO  3314  7S$Z1  213440  231434  774*0  *5014  ,*  <.«•*•)  * 

*03  3320  70*2.2  122123  141637  44405  91517  to 

•UA  3324  m?3  30302  90016  10402  10016 

903  1311  7Jfi£3  204104  223020  74444  *1900  *T0  t  «  W  “>  .  F  #0  t  -  ? 

SUA  1332  70024  53006  70346  14000  2942*  /Ol*« 

90S  3333  70024  71*00  02610  26166  10370  PC  L  ^ 

HUL  1334  70024  *4004  112000  11204  40600  1 « t  -  » 

*06  3  334  70024  11404  4  1  2152  5  42  0  4  4  4  5  325  **  6  jQu?  C 

•UA  1334  70024  141117  1*1422  11217  *0502  K> , 

90S  3140  75025  2092*  *0140  7*2*  14*00 

COO  134*  70025  1  44211  21341*  70431  77456  ’  /«,  /JJ  *  r." 

903  3353  T0*±6  102324  122015  17404  *4435  17, 31*'-'  *iC*>«*  ‘  * ’* 

COO  1154  70627  *5710  0*42*  25010  120*6  jl.JJi.rr  ?*TVA»*  AU 

•03  3364  70827  213039  233540  77414  0444*  jOC'S* 

HUl  1340  Totta  121452  141704  44142  51424  *4  L~la  TQ  C+ im 

HUL  3376  73529  05345  105054  32025  140*4  f  1 %  *  S’  ^  =  2  90  *  00 

NUL  1364  70030  02147  102101  30227  172*1  ^  llTl*  7  SOL*  7 

006  3360  70430  *00652  220605  72532  74**5  ,  JJ  *.  *  ■  -  >  £ Q  <-  *  fcr 

90S  3391  7044J  75357  45106  20437  154*0 

903  1402  7092.1  4212*  111014  31*0*  40714 

TOG  3407  70902  *5144  05053  24024  11093 

COO  34  1  3  7040  2  2114  5  7  2  1  3201  77*97  04  72  1  4  *1  J0"  *Q  ,  t  &  <-  ~  ~ 

HUL  3414  7592.3  101614  121940  1*474  4*1*0 

HUL  1425  70904  44*52  114929  19212  42129 

COO  1420  709115  203420  221139  740*0  01044 

•03  1412  70905  91605  111520  19405  4*520 

90S  1439  7 0925  230030  5741  9203*  1441 

90S  1442  ^  134115  151044  44249  1*12* 

HUL  1449  70907  4221*  *1420  1*71*  227*0 

006  1453  7 flaO 7  190250  205953  *0570  7*593 

NUL  3450  70900  44322  114027  1*002  42027 

HUL  144*  70409  41212  110422  *1112  401*2 

•UA  146*  70404  120021  140*31  41701  *0731  ll  L*  S~~  PO  l  •  **-  * 

ID  3444  70409  205*45  22*144  7*405  02424  4-  i.3  £  .  S’ 

COO  1472  70409  215141  14454  *4021  *344  0  »>T  /=~>  90L^U 

HUL  *473  70910  04350  10301*  Jl*10  1024*  j rr  *  \T'*L  5  U 

SUA  1473  70910  111710  11142  3  41030  400*3  *S  i  '  pb.iTO  17* 

TO*  147*  70910  202515  222020  71535  77100  U  ">>  >  *0 

BOS  1477  7Q9I0  213719  2222*9  77034  0O569  ,  '  ^ 

90S  5477  70410  232139  1*21  *6044  461  ^ 

ro*  1470  70911  2954  11051  1744  4731 

903  1400  70921  *1400  71104  10040  2*044 

•OS  1401  704J1  01*40  100711  31000  1*413  %  _ _ 

COO  1402  70411  112600  110314  41200  4*494  #  T>r  /  -  f./Vfc.r  ifO 
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ali ability  in  identifying  REPs  was  obtained. 

Fig.  C-3  shows  how  each  microfiche  is  layed  out  and  Table  C-l 
is  a  listing  of  the  individual  microfiche  by  Rev  Number,  Date,  UT 
of  Rev  in  hrs  and  secs,  and  notes  on  electron  and  proton  precipitation 
events.  Approximately  1212  microfiche  were  examined  three  times  by 
one  individual  (the  author)  and  once  by  another  graduate  student 
(T.  Kelley)  for  comparison.  A  list  of  358  outer  zone  transits  were 
identified  for  further  discrimination  by  pitch-angle.  313  outer  zone 
REP  revs  were  confirmed  by  the  pitch-angle  analysis  and  are  listed 
in  Table  C-2.  The  remaining  case  was  eliminated  by  the  loss-cone 
criteria  previously  described,  bad  data  drop  outs,  low  flux  values, 
and  spike  events.  This  table  also  subjectively  groups  the  revs  into 
173  'events'  by  combining  consecutive  revs  with  strong  precipitation. 
Since  each  revs  last  3  hrs  and  there  was  no  AE  data  set  available  to 
determine  substorm  duration  this  grouping  maybe  off  as  much  as  25%  in 
total  event  number. 
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Table  C-2.  Revs  with  strong  precipitation  >  235  keV  electrons 
by  EVENT  (one  or  more  consecutive  revs  with  strong 
precipitation),  DATE  (last  digit  of  year,  month, 
day),  Readout  REV  (first  letter  of  station,  number), 
LAT  (North,  degrees).  East  LONG  (degrees),  L-VAL, 

MLT  (hrs),  j(e/cm2  s  str  keV)  =  differential  flux 
value  of  precipitating  electrons  from  235  keV  to 
highest  energy  channel  at  which  strong  precipitation 

p 

occurred,  J(p/cm  s  str)  =  integral  flux  of  strong 
diffusion  precipitation  of  protons  per  energy 
channel  if  any,  80  p+  L-VAL  =  L  value  of  80  keV 
proton  precipitation  if  any,  strong  precip  12,  33, 
L-VAL  =  check  if  12,  33  keV  electron  channels  had 
strong  precip  and  L  value  of  12  keV  channel,  UT  (sec) 
time  at  which  235  keV  electron  precip  occurred. 
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Appendix  D  -  Computation  of  Electron  and  X-ray  Energy  Deposition 
Profiles 

The  precipitation  of  energetic  electrons  into  the  atmosphere 
is  simulated  by  a  finite  difference  numerical  solution  to  the  multi¬ 
dimensional  Fokker-Planck  diffusion  equation  developed  by  Malt  et  al . 
(1968).  The  technique  was  slightly  modified  to  handle  other  than 
isotropic  pitch  angle  distributions  by  Spjeldvik  (1974)  but  this 
modification  was  not  used  here  since  the  distribution  over  the  down¬ 
ward  hemisphere  at  300  km  is  isotropic  for  all  cases  considered. 

The  conversion  of  the  Malt  et  al .  (1968)  program  to  the  UCLA  computer 
by  Spjeldvik  (1974),  however,  was  used.  Spjeldvik  (1974)  contains 
the  computer  listing. 

The  Malt  program  solves  the  following  final  version  of  the 
Fokker-Planck  equation  as  derived  by  Chappell  (1968):  ox^/oc j  f -PA- 


where  Table  D-l  defines  the  symbolism.  Notice  that  the  stochastic 
averages  <  (AE)2>  ,  <(AaAE)>,  and  ^(AEAr)>  (small)  and  higher 
order  terms  have  been  eliminated.  Steady  state  (^-  *  0)  is  assumed 
since  precipitation  is  steady  on  times  scales  much  longer  than 
atmospheric  transit  times  for  precipitating  electrons.  Also,  radial 
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Table  0-1  Nomenclature  (Spjeldvlk.  1974) 


f  =  electron  distribution  function  averaged  over  cyclotron  phase 
t  =  time 

c  =  speed  of  light 
R^  =  radius  of  the  Earth 

L  =  L-value,  field  line  equatorial  distance  in  units  of 
E  =  electron  energy  in  rest  mass  units 
a  =  local  pitch  angle 
aQ  =  equatorial  pitch  angle 
r  =  electron  gyro  radius 
s  =  length  along  the  field  line 
q  =  electron  charge 
mg  *  electron  mass 

e  *  the  natural  logarithm  base,  2.7182... 
n.  =  number  density  of  atmospheric  species  i 
Zj  *  atom  number  of  atmospheric  species  i 
lj  *  average  excitation  potential  of  atmospheric  species  i 
*  minimum  scattering  angle  corresponding  to  maximum  impact 
parameter  (for  thermal  plasma,  the  Debye-length;  for  the 
neutral  atmosphere,  the  atomic  radius) 
oi  -  (e  (E+2)/^)172 

^  y  denotes  stochastic  ensemble  average. 
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diffusion  is  ignored  (£  =  constant). 

The  program  outputs  the  amount  of  average  energy  deposited 
per  unit  volume  at  each  altitude  below  the  incident  altitude.  Also, 
assuming  an  average  of  35  eV  is  required  to  produce  each  ion  pair, 
the  rate  of  ion  pair  production  is  output  versus  altitude. 

The  complete  altitude  profile  of  energy  deposition  by 
bremsstrahlung  is  computed  from  the  Berger,  Seltzer,  Maeda  (1974) 
Monte  Carlo  simulation  method.  This  work  is  essentially  an  extension 
of  the  earlier  (Berger,  Seltzer,  Maeda,  1970;  Berger  and  Seltzer, 
1972)  Monte  Carlo  simulation  methods  by  this  same  group.  Basically, 
a  computation  is  made  of  10,000  electron  trajectories  for  each 
incident  electron  energy  from  2  keV  to  2  MeV  for  electrons  injected 
into  the  atmosphere  at  300  km  altitude.  The  trajectories  are 
histories  of  individual  primary  and  secondary  electrons  as  they 
coulomb  scatter  (elastic  by  nuclei  collisions,  inelastic  by  orbital 
electron  collisions)  off  neutral  constituents  and  deposit  energy  at 
various  altitudes.  In  addition,  the  production  of  bremsstrahlung 
quanta  in  these  interactions  and  their  compton  scattering  and  con¬ 
sequent  photoelectron  absorption  of  the  bremsstrahlung  quanta  in  the 
atmosphere  is  followed.  The  Berger,  Seltzer,  Maeda  (1974)  results  of 
energy  deposition  as  a  function  of  atmospheric  depth  (mass  thickness) 
by  bremsstrahlung  from  monoenergetic  incident  electron  beams  for 
wide-area  precipitation  flux  isotropic  over  the  downward  hemisphere 
Is  given  In  Table  0-2  in  the  form  of  energy  deposited  per  unit  mass 
thickness  (ev/g  cm  )  per  incident  electron  energy.  From  this 
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TABLE  D- 


* 


.  Energy  deposition  function  for  Bremsstrahlung,  for  the 
case  of  uniform  wide-area  precipitation  of  an  electron 
flux  isotropic  over  the  downward  hemisphere.  The 
quantity  given  is  ABR  (Zm/T0),  in  units  of  cm2gnH. 
Number  in  parenthesis  indicates  powers  of  ten. 

(Berger,  Seltzer,  Maeda,  1974) 
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starting  point,  we  proceeded  to  work  towards  the  final  product,  a 
profile  of  ion  pair  production  by  bremsstrahlung  from  an  incident 
electron  energy  spectrum. 

The  CIRA  (1965)  reference  atmosphere.  Fig.  D-l ,  was  used  to 
scale  atmospheric  depth  to  a  representative  altitude.  The  same 
reference  atmosphere  was  used  to  compute  the  atmospheric  density 
(g/cm)  versus  altitude  (km).  At  each  altitude  the  density  was 

O 

multiplied  by  the  energy  deposition  function  Ag^[keV/g  cm  ],  to 
yield  the  stopping  ratio,  BR[keV/cm,  column],  normalized  to  an 
incident  electron  of  a  given  incident  energy  (To).  The  resulting 
column  stopping  ratios,  Br  (altitude),  are  listed  in  Table  D-3  and 
graphically  represented  in  Fig.  D-2.  Due  to  the  irregularity  of  the 
curves  and  diminishing  contribution  compared  to  electron  production 
above  50  km  in  Fig.  D-2,  only  value  below  50  km  were  considered.  For 
every  5  km  of  altitude  between  15  and  50  km  (inclusive),  a  linear 
(log-log)  interpolation  of  incident  electron  energy  (To)  versus 
column  stopping  ratio  of  Bremsstrahlung  photons  was  constructed. 

For  example,  at  15  km  altitude,  Br  =  10x  where  x  =  1.018342561  log 
To  (keV)  -  11.85638142  for  initial  energies  To  between  200-500  keV. 

A  computer  program  using  the  linear  (log-log)  interpolation 
curves  as  described  above  was  written  to  calculate  the  production 
of  ion  pairs  versus  altitude  for  an  input  electron  differential  energy 
spectrum  Isotropic  over  downward  hemisphere  at  300  km  by  solving  the 
equation 

Prod  (alt)  =  rwr  //  Br  •  Dif  (E)  •  dnd  E 
IP 
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TABLE  D- 


•  Br(keV/cm,  column),  stopping  ratio,  versus  altitude  for 
given  Incident  electron  energies. 
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(WX)  M 
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FIG  D-2.  Br(^,E)  versus  altitude. 
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where  EMI  =  average  minimum  ionizing  energy  taken  here  to  he  35  eV. 
Dif  (E)  =  differential  electron  flux  at  300  km 
[#/cm^  sec  str  keV] 
dn  =  solid  angle  ~  2tt  here 
dE  -  1  keV  intervals 

Prod  (alt)  =  ion  pair  production  rate  at  each  altitude 
[#/cm^  sec] 

Also  output  is  the  total  energy  (keV)  deposited  at  each 
altitude  per  (cm  str  sec).  The  computer  listing  follows: 
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TABLE  0-4.  Bremsstrahlung  deposition  computer  listing. 
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